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INTRODUCTION 

This series of notes tries to lay the foundations for the interpretation of magnetic fields and 
fluxes, often in terms of equilibrium plasma parameters. The title, 'magnetic diagnostics', is 
taken to mean those diagnostics which are used to measure magnetic fields and fluxes using 
induction, or pick-up, coils. More specifically, what is often inferred is a question: "How much 
can we tell about a plasma given certain measurements of magnetic fields, and fluxes, outside 
that plasma?" I don't consider here diagnostics which measure the plasma current density 
distribution utilizing phenomena such as the motional Stark effect, or Faraday rotation; these are 
found in a series of notes on Plasma Diagnostics.. 

The measurements themselves are in principle simple, although in practice they are always 
complicated by unwanted field components, for example from misaligned pick-up coils. There is 
also the problem of allowing for image currents flowing in nearby conductors; dealing with these 
image currents becomes a large part of the problem. Including the effects of an iron core also 
leads to complications. 

Many people think the topic under consideration is boring, in that there is nothing new to do. 
You have only to read current issues of plasma physics journals to recognize that there is still 
much interest in the topic. For example, equilibrium and its determination, axisymmetric stability 
and disruptions are all of current interest, and all involve 'magnetic diagnostics'. The subject 
does appear to be difficult (students starting in the topic have a hard time). 

The layout of the notes is as given in the list of contents. Generally I have included topics which 
I have found useful in trying to understand tokamaks. Some basic concepts (inductances, fluxes, 
etc.) are included, because they are made use of throughout the notes. There is also a section on 
plasma equilibrium, in which the large aspect ratio, circular tokamak is described. The fluxes 
and fields from this model are used as examples for application of certain ideas in the remainder 
of the text. 

References I find useful include: 

o B. J. Braams, The interpretation of tokamak diagnostics: status and prospects, IPP 

Garching report IPP 5/2, 1985. 
o L. E. Zakharov and V. D. Shafranov, Equilibrium of current carrying plasmas in toroidal 
configurations, in Reviews of Plasma Physics volume 11, edited by M. A. 
Leontovich, Consultants Bureau, New York (1986). 
o V. S. Mukhovatov and V. D. Shafranov, Nucl. Fusion 11 (1971) 605. 
o V. D. Shafranov, Plasma Physics 13 (1971) 757. 

o L. E. Zakharov and V. D. Shafranov, Sov. Phys. Tech. Phys. 18 (1973) 151 
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o J. A. Wesson, in Tokamaks, Oxford Science Publications, Clarendon press, Oxford, 1987. 
o P. Shkarofsky, Evaluation of multipole moments over the current density in a tokamak 
with magnetic probes, Phys. Fluids 25 (1982) 89. 
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1. SOME CONCEPTS AND DEFINITIONS 



Maxwell's equations 



We are going to make extensive use of Maxwell's equations. In vector form, these are 



VxB = 




1.1 



VB = 



1.2 



3B 

VxE = 

at 



1.3 



VD = p 



1.4 



If charge is conserved we can add to these the continuity equation. We shall ignore the 
displacement current SD/dt, and take (I = (i , the free space value, inside a plasma. Without the 
last term Equation 1.1 is Amperes law. We have effectively restricted ourselves to assuming nj = 
n e = n, the charge neutral assumption, and that any waves have frequencies much less than the 
electron plasma frequency, with characteristic lengths much greater than the Debye length . We 
have not said E or V-E = 0. When considering plasma equilibrium we shall also assume the 
electron mass m e approaches 0. This allows electrons to have an infinitely fast response time. 

Pick-up or Induction Coils 

This is the heart of the matter. Magnetic fields are usually measured with pick-up or induction 
coil circuits. Changing the magnetic flux in a circuit generates a current; the direction of this 
current is in a direction such as to set up a magnetic flux opposing the change. The 
electromotence or voltage £ (V £ = -E, the electric field intensity) in Volts induced in a circuit 
equals the rate of change of flux N = pB ■ ndS in Webers per second, i.e. 



The flux can be changed either by changing its strength, changing the shape of the circuit, or 
moving the circuit. Then 



dN 



1.5 



£ = 



dt 




1.6 
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for any path 1, with n the normal to a two sided surface S. Applying Stokes theorem 

(f A • dl =Jn • V x AdS for any vector A) to the left hand side of Equation 1.6 gives Equation 

s 

1.3. Figure 1.1 shows the geometry of a coil used in applying Equation 1.6, Faraday's Law. The 
output signal must be time integrated to obtain the required flux. By taking a small enough coil 
the local field B can then be determined. This becomes difficult if very small scale variations in 
field exist, because the pick-up coils must then be very small themselves. The surface S includes 
any area between the leads; this is minimized by twisting them together. A hand drill is 
particularly useful for this. 



Coil 




voltage O ' Surface S 

produced f V- E_E^ 

across k 
leads \ 

leads ~ N Contour 1 

Figure 1.1. The contour 1 and surface S of a pick-up coil. 



Integration 



The time integration required to obtain the magnetic field B from the pick-up coil output e can be 
performed either digitally or by an analog circuit. 




'in 



8 

out 



Figure 1.2. A passive "QC" integration circuit. 

The simplest thing to do is to use a capacitor (C) and resistor (Q.) network, as shown in Figure 
1.2. The output voltage is given by 

de, 



out i O out _ in 



dt 



1.7 



with x = QC called the integrator time constant. The solution to this equation is 
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e out — e 



---Xj(e-*)e in (.'A 



For example, suppose at t = we start an input voltage £j n = £j n o sin(cot), so that the required 
integral is £j nt = £j n o (l-cos(cot))/co. The output from the passive circuit is (obtained using 
Laplace transforms) 



( ^ 

cox sin(cot) - cox cos(cot) 

^out — ^inO \7 7 T\ 

e^(l + (cox) 2 ) (l + (cox) ) 



1.9 



Now consider two extremes. First, if cox » 1 and t « x we have 

e out =-^(l-cos(cot)) 1.10 
cox 

That is, £ ou t = 1/t times the required integral. In this limit we have integrated the input signal. If 
cox « 1 and t » x, then £ out = £i n . 

As an example, we show in Figure 1.3 the output from the passive integrator ("integrator 
output", dotted line) for a sinusoidal voltage input of IV at a frequency of 100 Hz ("coil input", 
solid line), with an integrator time x = 0.1s. The exact integral ("field") divided by x is shown as 
the broken line. The integral is only performed accurately for times t « x; as the pulse proceeds 
there is a "droop", and significant errors result. We can imagine the curve "field" represents a 
specified magnetic field time history B = Bo(l-cos(cot))/(cox), with Bo = x/(nA) T, and B/x is 

plotted. The curve "coil output" represents the unintegrated output from a magnetic pick-up coil 
with area nA m 2 , (n turns each of area A), which becomes the input voltage to a passive 
integrator £j n = sin(cot). Finally the curve "integrator output" represents the output from the 

passive integrator, which we would interpret as the original magnetic field. 

A common situation is that the required signal from the pick-up coil has a low frequency 
component of angular frequency C0q, and superimposed upon this is a higher frequency unwanted 
"noise" signal of angular frequency C0i. By carefully choosing the time constant x of our passive 
integrator so that coot « 1 (£out = £in) but COiX » 1 (integration) we filter the noise, leaving the 
required slowly time varying voltage. As an example, Figure 1.4 shows the passive integrator 
output £ ou t (dashed line) for an input voltage £i n (solid line) comprising a slow (coo = 10 rs _1 , £i n o 
= 1 V) and fast (C0i = 2xl0 3 rs _1 , £i n i = 0.2 V) component. The time constant x = 0.01 s, so that 
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(Qq% = 0.1 (« 1) and a>ix = 100 (» 1). The output voltage is filtered, as required. The dashed 
line shows the exact integral divided by x, for comparison. 

signal 




Figure 1.3. The input ("coil input", solid line) sinusoidal voltage with f = 100 Hz and 
output ("integrator output", dotted line) of a passive integrator circuit with x = Q.C = 
0. 1 s. The (exact integral)/x is denoted by "field", the broken line. 




time (s) 

Figure 1.4. The output of a passive integrator circuit used as a filter. An input voltage 
(solid line) with summed sinusoidal voltages is smoothed to give the dash-dot line. 
The exact integral divided by x = Q.C is shown as the dashed line. 

A more common system to perform the time integration is an active integrator, but in many cases 
an input filter consisting of a passive integrator is still used. Active integration is performed 

using a circuit such as shown in Figure 1.5; the output voltage e out = [ £ in dt . The example 

h 

shown grounds one side of the coil. A useful feature shown is the integrator gate, which defines 
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the time ti the integration starts. On tokamaks this gate is often used to help reduce errors from 
misaligned pick-up coils. For example, tokamaks have a large toroidal field and a much smaller 
poloidal field. Therefore if the pick-up coil used to measure the poloidal field is misaligned even 
by a small amount, the resulting component of the toroidal field which is picked up (as dB/dt) 
can be significant. However the toroidal field usually evolves on a much slower time scale than 
the poloidal field, and in fact it is usually time independent at the time the poloidal field is 
initiated. Therefore the integrator gate can be opened when the toroidal field is time 
independent, and therefore the induced voltage in the misaligned pick-up coil is independent of 
the toroidal field. 



If the data is digitized, integration can be performed numerically. Sufficiently fast systems now 
exist for "real time" integration; the integration can be performed in (is so that integrated signals 
suitable for real time feedback control can be obtained. A p bit digitizer has a resolution of 1 part 
in 10P, e.g. an 8 bit system has a resolution of 1 in 256, while a 10 bit system has a resolution of 
1 part in 1024. This can be a limitation if we intend to investigate large but low frequency 
magnetic fields in the presence of small, high frequency fields. An example is that of trying to 
measure the equilibrium poloidal field in the presence of Mirnov oscillations. The pick-up coil 
output is dominated by the voltage produced by the time derivative of the small but high 
frequency component. Avoiding saturating the input by the higher voltage, high frequency 
component means that the resolution of the low frequency fields is now restricted. If we want to 
use the full capability of the digitizer in recording the lower frequency fields, then the solution is 
to filter the signal and only allow frequencies below a certain value to be recorded, i.e. use the 
filter described above with reference to Figures 1.2 and 1.4. 

Intuition suggests that if a time varying wave form is sampled sufficiently fast then the original 
wave form can be recovered. However, we must determine how close the samples must be, and 
how to interpolate between adjacent points. The sampling theorem provides answers to these 
questions. An original signal x(t) can be recovered from sample values x(nt s ), with t s the sample 
time, by locating sine functions at nt s with amplitudes x(nt s ). The signal x(t) can only be 




Figure 1.5. An active integrator circuit. 
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recovered if the signal bandwidth b < f s /2, with f s the sampling frequency = l/t s . If this is not 
done, aliasing occurs. 

If b > f s /2 then the high frequency signal can appear as a low frequency signal. The fact that 

spoked wheels in films sometimes appear to rotate backwards is a manifestation of aliasing. 
Aliasing can be avoided using a passive filter to remove the high frequencies f > f s /2. For 

example, sampling at 5 kHz (i.e. a sample every 0.2 ms) then an "anti aliasing" filter with x = 0.5 
ms can be used. 

Vector potential 

In describing plasma equilibria we shall make use of the vector potential A. It is related to the 
poloidal flux, and used to determine self and mutual inductances. It is defined through the 
equation 

VxA = B 1.11 
In cylindrical geometry (R,(|),z), which we shall use a lot, this is 

R R d(/) dz 
dA R dA 

_l 3(*A,) 13A P 
z R dR Rd</> 

Then the electric intensity E is proportional to the vector potential A whose change produces it: 

E = -^ 1.13 
dt 

From Equation 1.1 (ignoring D and p) we then have that A is given in terms of the current 
density j (the current per unit area) by Poisson's equation: 

V 2 A = -uj 1.14 

which has a solution 

A = iLf jdV =iL ^Idl 
J V r 



47! Jv r 471 r 
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where the total current I flows inside the volume V, the line element dl is along the direction of 
the total current I, and r is the distance from the line element to the point of interest. A useful 
example for us is the vector potential of a circular filament. This is used to represent windings 
(vertical field, shaping, ohmic heating) on the tokamak, and elements of the plasma current itself. 




Figure 1.6. The geometry used to evaluate the vector potential of a circular filament. 

Consider a circular filament of radius Ro, with current I in the (]) direction. must be 
independent of (]), so choose the point of interest P in the (X,z) plane of Figure 1.6, where (]) = 0. 
Pairing equidistant elements dl shown in thickened lines at +(]) we see the resultant is normal to 
(R,z). Therefore only consider the component dl^ of dl in the direction normal to the plane (R,z); 
dl^ = Rocos((]))d(|). The radial distance r from the point P to the element is given by 
r 2 = z 2 + +R 2 -2R Rcos(<f>). Then 



1.16 



Ak r r 2^J^2 + jR 2 +z 2_ 2 ^ /?cos( ^ ) y4 
Far from the loop (i.e. a small loop) we have rg = (R 2 +z 2 ) l/2 » Rq, and the integral becomes 



_ fill f /^cos(^) 

ln 'o r o V 
RRlluI _ (Mxr) 
Ar Am- 3 



RRo cos {(f>) 



d(j) 



1.17 



Here we have written the magnetic moment of the loop M= 7rR Q I , directed upwards. 



If the loop is not small, then let § = n + 20, so d(|) = 2 d0 and cos((|)) = 2 sin 2 (0) - l, and we obtain 
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MM 
n 



"A 

\ 



(2sin 2 (6')-iy^ 



((# + Rf + z 2 - 4R Q R sin 2 (6)f 



1.18 



This can be re- written in terms of K(k 2 ) and E(k 2 ), the complete elliptic integrals of the first and 

dd 



2. £ 
2 I 2 

second kind, (E(m) = J(l-msin '(0))d0, K(m)=\ . 

o ^y(l-msin 2 (0)) 



as 



_Ml(Ro) y > 



£ 2 ^ 

1 -— ^ 

2 J 



£ 2 = 4fl /?[(/? + fl) 2 +z 2 ]' 



1.19 




Figure 1.7. Cylindrical (R,(|),z) and quasi cylindrical (p,(0,(|)) coordinate systems. 
Both are right handed 

Going into a quasi-cylindrical coordinate system (p,co,<f)) shown in Figure 1.7, based on the 
current path, then 



R =R - pcos((0) 
z = psin(o)) 



1.20 



Expanding k we find that keeping terms of order p/Rg then k 2 ~ 1 . We must go to higher order 
(to find k < 1) because K(l) = <». Keeping terms of order (p/Rg) 4 we find 

2 1 p 2 p 3 cos (a) p 4 (l + 2cos(2tf;)) 
k ~1 j ) j . Using the expansions for E and K found in the 



ARl Rl 4 



K 



16 



Handbook of Mathematical Tables, Abramowitz, Dover Publications, and after some 



rearranging, we obtain up to order (p/Rg) 2 E 



( 


r 2R )) 


1 + ln 


V 


v p J) 



and 
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K~ In 



v P 



P r \ P 

-cos(ft>) + - 



2R, 



f 


r 2*: 




In 




+ 


V 


v p J 





+ — (l + 4cos(2<»)) . Finally we can write an 

2 J 



expression for the vector potential near the loop keeping terms of order (p/Ro): 



Mo 1 



2k 



In 



Ml 
V P 



-2 



+ 



p cos(cd) 
\ 2 



( f 
In 



v P 



1.21 



Figures 1.8 and 7.9 show a comparison between the approximate (Equation 1.21) and exact 
(Equation 1.19) solutions for a case with current I = l/|i, and a major radius Ro = 1.0 m. For the 

approximate solutions we show results both for zero order (i.e. neglecting terms proportional to 
p/Ro) and including first order corrections (i.e. including terms proportional to p/Ro)- Figure 1.8 
shows a cut in the plane of the coil (z = 0), while Figure 1.9 shows contours for = 0.15, 0.2, 
0.3 and 0.4. We see that it is sufficient to consider only the zero order terms, that is the terms 
proportional to p/Ro can be neglected. In fact in Figure 1.8 the effects of including the first order 



terms cannot be detected. The simple approximation ~ 



ln^ 

V V p 



is excellent. 



J 



The field components are given, in the original circular coordinate system (R,(]),z), as 



R dz 2k 



1 



2 2 

+ z 



b i^K) pj i r Ri-^-e E ~\ 

where we have made use of 

dK _ E _K.dE_E_E. 
dk~k(l-k 2 ) k'dk'k k 9 

dk _zk 3 dk k k 3 k 3 

dz ~ 



4R R'dR 2R 4R 47^ 
On the axis R = we have Br = and B = ■ 



iR 2 Q +z 2 f 



1.22 



1.23 
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0.5 0.8 1 1.2 1.4 

major radius R (m) major radius R (m) 

Figure 1.8. The exact (solid line), the approximate Figure 1.9. Contours of vector potential 
zero order (dashed line) and the approximate first from the exact solution (thin dark line), 
order (dotted line) solutions in the z = plane for the zero order approximation (thick gray 
the vector potential from a circular current filament line) and including first order corrections 
with current I = 1/(1. (intermediate thickness dark line) for a 

circular current filament with current I = 
1/(1. Contours of potential equal to 0.15, 
0.2, 0.3 and 0.4 are shown. 



Mutual inductance 



We shall use mutual and self inductances, often between circular filaments. They can be used to 
derive the vertical field necessary to maintain a plasma equilibrium, and to analyze axisymmetric 
instabilities. We are interested in the relationships between mutual and self inductances and 
fluxes, and how to write energy in terms of mutual and self inductances. 

The mutual inductance M12 between two circuits is defined as the flux N12 through circuit 1 
produced by unit current in circuit 2. Then 

M n = ^dl x 1.24 
with A2 the total vector potential due to current 2 in circuit 2. That is, 
4ff Jl i Jl 2 r 

The electromotance through circuit 1 due to a current I2 in circuit 2 is 
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e =M 4k 126 

1 12 dt 

The total energy (in a volume V) associated with two circuits is 

1 1 1 121 

= — f B?dV+ — f bUv + —\ B, »B 2 dV 

The first two terms represent the energy required to establish the currents Ii and I2 producing the 
fields Bi and B2 in circuits 1 and 2. The third term is the energy used in bringing the two circuits 
together. This mutual energy between the two circuits W12 is 

W u = M n IJ 2 =-lB l B 2 dV 1.28 

Self inductance 

The magnetic energy density W of a single circuit carrying current Ii is used to define the self 
inductance I4 1 of the circuit: 

W = —[B 2 dV = -L.lf 1.29 
2ju Sv 2 M1 1 

To maintain the current Ii a power source must, in each second, do an amount of work 
dt 

(because £ = dN/dt) in addition to working against resistance Q. The stored energy per second in 
the magnetic field equals dW/dt, so that 

/*= I , 1 ,i 1.31 

1 dt 11 1 dt 

i.e. 

N n = L n I x 1.32 

Therefore we can also define the self inductance of a circuit through the change in flux linking 
that circuit when the current changes by one unit: 
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dN 1L= dL 
1 dt 11 dt 

Poloidal flux 

Suppose that a system consists only of toroidally wound loops producing only poloidal fields. In 
a cylindrical coordinate system R,(|),z shown in Figure 1.7 is also the 'toroidal' angle in a quasi 
cylindrical coordinate system) nothing depends on the angle (]). Then 

M n = 2xR l ^ 1.34 

^2 

and A has only a toroidal component A^. In this case the fields are given by (B =VxA): 

B = 

dz 

B^=0 1.35 

B _ i a K) 

z a/? 

These poloidal fields are also expressed in terms of the transverse (poloidal) flux function 
x ¥: ^(R^z) = constant defines the form of the equilibrium magnetic surfaces, proved later: 

B = — (V^xe.) 1.36 

with e§ a unit vector in the toroidal ((])) direction, so that 
1 d*¥ 

1.37 



B =- 

z 2kR dR 



R 2M dz 

But we know that we can write, in terms of the vector potential for our toroidally symmetric 
system (d/d§ = 0), 

B R = 



dz 

, . 1.38 

i afc*,) 



B 

z R dR 



That is, the poloidal flux can be written as (with subscripts implied but not given): 
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¥ = 2xRA d = MI= 2tu\ R BMR 

<p Jo z 



1.39 



That is, for the system we are considering, the poloidal flux at a position R is simply the vertical 
field B z integrated across a circle of radius R. Note that sometimes in the literature the flux 

function \\f = *P/(2tc) is used. 

As an example we show in Figure 1.10 the poloidal flux \\f produced by a single circular 
filament (see the section on vector potentials for the derivation of A^) of radius Ro = 1 m, current 
I = 1/jxO- Results are shown in the plane of the coil (z = 0). The exact results are shown as the 
solid line. Also shown are two approximate solutions; the very near field solution and the far 
field solution. Near the current ("near field") we can write 



y=RA,~ 



In 



( f 
In 



8 ft 



v v p 



p \ cos(<w) 



f f 

In 

v v p 



1 



1.40 



The very near field solution is the zero order in p/Ro term, i.e. y/ = 

from the loop ("far field") we can write 
MR 2 



In 



( 






In 


-2 


V 


v p J 





(R 2 +z 2 f 



1.41 



Far 



We see that neither the very near or far field solutions are good for distances of about half a coil 
radius from the coil itself. However, if we use the full expansion expression including terms of 
order p/Rg (Equation 1.40) then the results are very near to the exact solution. This is shown in 

Figure 1.11. 




0.5 



R(m) 



1 .5 
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Figure 1.10. The poloidal flux \\f in the plane z = for a circular current loop, 
radius 1 m, current I = 1/\Xq. The solid line is the exact solution, the long dash line 
is the far field solution, and the short dash line is the very near field solution (zero 
order in p/R terms only). 
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0.5 0.8 1 1.2 1.4 

Radius R (m) 

Figure 1.11 The poloidal flux \\r in the plane z = for a circular current loop, 
radius 1 m, current I = I/jiq- The solid line is the exact solution, the short dash line 
is the first order expansion in p/Ro solution keeping terms of order p/Ro, and the 
long dash line is the very near field solution (zero order in p/R terms only). 

We can also compare the zero order, first order and exact solutions by plotting contours of \\f in 
(R,z) space. This is done in Figure 1.12, for the conditions described in the caption of Figure 
1.11. We see that it is very important to include the first order in p/Ro terms; even then the 

solution contours are significantly different from the exact solution contours. This is very 
different from the case of the vector potential, where the zero order solution was accurate. 
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0.6 0.8 1 1.2 1.4 

radius R (m) 

Figure 1.12 Contours of flux \|/ from the exact solution (thin dark line), the zero 
order approximation (thick gray line) and including first order corrections 
(intermediate thickness dark line) for a circular current filament with current I = 
1/(1. Contours of flux equal to 0.15, 0.2, 0.3 and 0.4 are shown. 

Field lines and flux surfaces 



Before starting on TOKAMAK EQUILIBRIUM, I want to discuss the field line equation, and 
apply it to a case very similar to a tokamak. The field line equation is 



dx dy dz 
B ~ B ~ B 

x y z 



or, in cylindrical geometry (the system (R,(|),z) of Figure 1.7) 
dR = dz = RB £ 

B R B z B & 



This is easy to see in two dimensions, as illustrated in Figure 1.13. 



1.42 



1.43 
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Z 
A 




>X 



oc= dz/dx = 

Figure 1.13. A field line made up from two components B x and B z . 

If the length along a field line is 1, then these equations are equivalent to dl/B = constant (B is the 
magnitude of B). Now a magnetic surface is defined by an equation V|/(R) = constant; we use \|/ 
because it will turn out to be the poloidal flux that is constant on a surface. The condition that all 
lines of magnetic force lie upon that surface V|/(R) is written as 



V^«B = 



1.44 



because V\|/ is the normal to the surface, so the equation says there is no component of B 
perpendicular to the surface. Now we describe the fields B through the vector potential A, so 
that 



dz 



R R d0 



„ dA R dA. 
* dz dR 

z R dR Rd(/) 



1.45 



For axial symmetry, nothing depends on (]) (we will also call this toroidal symmetry). Now take 
y(R,z) = RA, 1.46 



Using Equation 1 .45 in 1 .46 gives 

d(RA,) [i dA 
R component: — 1 — — L — — ^ i = - 

dR IR d(/> dz J 

i aK) f3A fi at] n 

^component:- — L —-^ r 



R— A +A ' 



dR 



>) dz 
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a(*A,Yia(*A,) iaA s l a^r ^ 

z component : — — — - = — ^ R—^ + A, i 1 .47 

F a z [/? a/? a^ J aA a/? ') 

where the RHS of each equation represents the result for dld§ = 0. Therefore, with the assumed 

symmetry, the § component is zero and the R component and z component cancel. Therefore our 
assumed form for the surface (Equation 1.46, y/(R, z) = RA^) ensures that V^«B = 0, i.e. it 

ensures that all field lines lie on that surface where \\f = RA^ = constant. In the case of toroidal 

symmetry the z component of Equation 1 .45 gives 



_l8K)_ 13(0 
z R dR R dR 

R j R 

i.e. y/ = \ B z RdR = — \ 27iB z RdR 1 .48 

o 2;r o 

l ¥? v 

u BzRdm 



2^ J 2^: 

where the total poloidal flux *P is the integral of the vertical field B z through the circle we are 
considering. 

An example 

I want to consider what the magnetic surfaces look like starting with a single filament in the § 
direction (a 'toroidal' current), then adding a uniform 'vertical' field in the z direction, and 
finally adding a filament in the z direction to produce a 'toroidal' field. This is meant to 
approximate a tokamak equilibrium, but note we are not specifying that the conductors are in 
equilibrium yet. 

Consider a circular filament, as shown in Figure 1.6, with current I^. The field lines lie in a 

surface (the flux surface) defined by \\r = constant, and we have shown previously that to zero 
order in the normalized distance p/Rg from the filament 



( r n „\ \ 



In 



M»,_ 2 

V p J J 



= constant 1 .49 



lit 
where 

p 2 =(R Q -R) 2 +z 2 1.50 
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i.e. the magnetic surfaces are described by circles with p = constant. 

Far away from the loop we have shown that the field looks like that due to a dipole with moment 
M: 

M = tcR 2 V 1.51 
Then the equation for the magnetic surfaces becomes 
MRl 

w = v = constant 1.52 

(R 2 + z 2 f 

Now we add a uniform field B z o in the z direction. This has a vector potential given by 
A„=2f 1.53 

The magnetic surfaces are now given by 

r(a^ + A^ )= constant 1.54 

Finally we add a filament up the z axis, which produces a field 1/R in the (]) direction (a 
'toroidal' field. The vector potential due to this filament is 

^ _ _JhL. \n(R) 1.55 
Itc 

Because this has only a z component, it does not affect the result (Equation 1.54). The results are 
shown in Figure 1.14 and Figure 1.15 for a vertical field (the form used is relevant to that 
required to maintain a circular tokamak in equilibrium) 

z 4xR V V a J 2) 

with A = 2, Ro = 1 m, and a = 0.2 m. In the figures the exact form for A^ from the circular 
current filament (written in terms of elliptic integrals) is actually used. The positive current 

produces a field downwards (in the -z direction) at the outer equator (z = 0, R > 1 m). Adding a 
positive vertical field B z cancels this field at some point, producing a point where B z (z=0) = 0. 

This is called an X point; the flux surface through this point is called the separatrix. Note that, 
with the negative uniform vertical field applied, there is no inner X point. However, if a vertical 
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field oc 1/R was applied the more negative B z at the inner equator would cancel the positive B z 
from the filament, and an X point would appear 




0.6 0.8 1 1.2 1.4 0.6 0.8 1 1.2 1.4 

major radius R (m) major radius R (m) 

Figure 1.14. Contours of constant flux \|/, Figure 1.15. Contours of constant flux V|/, 
with positive B z . with negative B z . 

x 




Figure 1.16. The geometry used for calculating the flux surfaces of a straight filament in a 
quadrupole field, z is into the plane of the figure. 
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For a second example, consider the field lines and flux surfaces resulting from a quadrupole field 
applied to a single filamentary current, in straight geometry. In a straight system the field lines 
will lie in a surface defined by constant vector potential A (i.e. \|/ constant, with R => <*>). The 
single ("plasma") filament of current I zp lies at the origin of a rectalinear coordinate system 

(x,y,z) shown in Figure 1.16, where z will approximate the toroidal direction in a toroidal 

system. The vector potential at P is given by Equation 1.55 as A = - — zp ln(r), with r the 

Ik 

radius from the plasma filament to a point P. There are then four additional filaments, at a 
distance d from the plasma filament, with currents I zq alternatingly + (into the plane) and - (out 

of the plane). The i th additional filament (i = 1 to 4) has a vector potential A . = - — <? ln(r ; ) in 

Ik 

its own local coordinate system. Transforming to the coordinate system (r,0) we obtain 



The resulting contours of constant A z are shown in Figures 1.17a through 1.1 7d for the ratio 
Izq/Izp = 1-0 to 2.0. In the examples d = 0.25 m. The equi-distant contours shown are different 
from frame to frame. Clearly approximately elliptic cross sections are obtained; in tokamaks 
elliptic surfaces are produced by applying quadrupole fields. A separatrix appears and defines 
the last closed vector potential surface (a straight circular system is the only one without a 
separatrix). 

We can derive an analytic expression for the shape of a given flux surface. By expanding 
Equation 1.57 to third order in r/d, we find 



With I zq = the surfaces are described by r = constant = a, i.e. a circle. If we assume the surfaces 
with I zq 4- are described by 




ln(d 2 + r 2 -2drsin(0)j' 2 j 
ln(d 2 + r 2 + 2Jrsin(0)) 1/2 J 



1.57 




1.58 



r = a + A 2 cos(2 0) 



1.59 
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Figure 1.17a. Vector potential contours with Figure 1.17a. Vector potential contours with 
Izq/Izp = 1- The quadrupole filaments are at I Z q/Izp = 2. The quadrupole filaments are at 
± 0.25 m. ± 0.25 m. 

then Equation 1.58 gives, after expanding in A/a, 
2k W 2k 



-A, = 



K 



+ r d ) k-} 



K 



1.60 



For a/d « 1, A/a « 1 we can ignore the last term (°c cos 2 (20)). Then we ensure that 

A 7 = const = - — zp ln(a) by setting the coefficient in front of the cos(20) to zero, i.e. 
2k 



* = 2 

a 



(l > 



1.61 



The elongation of the surface is then 
1 + 2 



height 
width 



f a 



1-2 



M. 



f a 



1 + 4 



<1 > 



1.62 
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Figure 1.18a and 18b show the computed distortion to a particular surface (A = constant) for d 
= 0.25, I zq /I Z p = and 2.0. With I Z q/I zp =1.0 the value of height to width is measured to be 1.28, 

as compared to the value of 1.25 derived from Equation 1.62. 
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Figure 1.18a. Contours of vector potential Figure 1.18b Contours of vector potential 

with I Z q/I zp = 0. The larger contour is taken with I Z q/I zp = 2. The contours have the same 

as a reference in determining the distortion flux values as those shown in Figure 1.18a. 
produced by an applied quadrupole field. 

Circuit equations 

For some applications we will consider the plasma as a lumped series resistance and inductance, 
coupled to other circuits (including a conducting vacuum vessel) by mutual inductances. Figure 
1.19 shows how this is represented. 

The equation for circuit 1 consisting of a series self inductance Ly and resistance coupled by 
mutual inductances My to other circuits i, is 



1.63 



The sum over the mutual inductances is for i ^ 1 because Mn = Ln is brought out separately. If 
the circuit is closed (short circuited), then 8i = 0. If the circuit is open, or connected to a high 
input impedance, then h = and £i = d/dt(EMi iIi). The plasma is sometimes represented as one 
series inductance-resistance circuit, or sometimes as a number of such circuits in parallel, all 
short circuited together. The vacuum vessel is similarly represented as a number of paralleled 
resistor-inductor circuits, which can be open circuit (a vessel with an insulating gap) or short 
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circuited (no insulating gap). For a perfectly conducting plasma the plasma series resistance 
elements are set to zero, so that the flux enclosed by the plasma loop is conserved. These models 
are very useful in analyzing axisymmetric stability problems. 




vessel/shell (switch open = gap) 
Figure 1.19. Representation of circuits coupling to a plasma. 
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2. SOME NON STANDARD MEASUREMENT TECHNIQUES 

Hall Probe 

Other techniques (than pick-up coils) are used to measure magnetic fields. The most common 
alternative is a Hall probe, shown in Figure 2.1. A semiconductor is placed in a field B, and a 
current I driven perpendicular to B. The current carriers experience a Lorentz force, producing a 
charge buildup in the direction perpendicular to both B and I. The resulting charge buildup 
produces an electric field which cancels the magnetic force. This electric field is measured by 
electrodes. This was iscovered in 1879 in Johns Hopkin University. 

Drive j. Charge build 
up perpendicular to B 
and j. Measure E. 



Figure2.1a. A Hall probe. 




Figure 2.1b. A Hall probe in use 
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Assume the electrons move inside a flat conductive strip in a magnetic field. Then the output 
voltage is: 



Where i is current, h is efficiency which depends on geometry, temperature, area. Theoretically 
the overall efficiency depends on the Hall coefficient, the transverse electric potential gradient 
per unit B field per unit current density. 



Specific problems include: susceptibility to mechanical stress, and temperature (of resistors). 

Faraday Effect 

It has also been proposed to use the magneto-optic effect (the Faraday effect) in fused silica 
single mode optical fibers to measure magnetic fields, and the electro-optic (Kerr) effect to 
measure electric fields. The Faraday effect is the consequence of circular birefringence caused 
by a longitudinal magnetic field. Circular birefringence causes a rotation F of the plane of 
linearly polarized light, given by 

F=V c j>H»d\ 2.2 



around a contour 1. No time integration is required. The Verdet constant V c ~ 5xl0~ 6 radA -1 for 
silica. Thus the rotation must be now measured. Another approach is to coat a fiber with 
magnetostrictive material and measure the strain effects, with the fiber as one arm of a Mach 
Zender interferometer. 



V H = hiBsin(a) 



2.1 




Figure 2.1c. A Hall probe at and angle to the field 



The Compass. 



Chinese 2634 BC, magnetite suspended on silk. 
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Flux gates 

This is intended for weak fields. A B-H curve below is shown below. 



Alan Wootton 




.\ Mj» kiil' (4cW Slmnjftt 



Figure 2.3 a A simple flux gate and the B-H curve 

An applied field H to the core induces a magnetic flux B = mH. For high B the material saturates 
and (I is very small. There is hysterisis, and the path is different for increasing and decreasing H. 
When the core is not saturated the core acts as a low impedance path to lines of magnetic flux in 
the surrounding space. When the core is saturated the magnetic field lines are no more affected 
by the core. Each time the core passes from saturated to unsaturated and backwards, there is a 
change to the magnetic field lines. A pickup coil around the core will generate a spike. Flux 
lines drawn out of core implies positive spike, lines drawn into core, a negative spike. The 
amplitude of the spike is proportional to the intensity of the flux vector parallel to the sensing 
coil. The pulse polarity gives the field direction. 

The core must be driven in and out of saturation by a second coil. The excitation current will 
induce a corresponding current in the sensor coil, but this can be allowed for. 
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Figure 2.3b. A driven flux gate 

A better approach is to position the excitation coil so that it will excite without affecting the 
sensor coil. i.e. excite the flux at right angles to the axis of the sensor coil. One can use a 
toroidal core with a drive winding and a cylindrical sensor coil.. 





Figure 2.3c: A flux gate with toroidal core with a drive winding and a cylindrical sensor coil 
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3. GENERAL FIELD CHARACTERIZATION 

Fourier components 



Suppose we want to characterize the tangential (subscript x) and normal (subscript n) fields on a 
circular contour of radius a.\. It is often convenient to express the results as a Fourier series: for 

the poloidal (0) and radial (p) fields outside a current I we can write 



i r 

B (0 =B T = -j^- ^1 + X K cos(n«) + 8 n mv{n(0) 
B p = B n = W K n cos(n«) + ju n sin(nO))] 



3.1 



3.2 



We are working in a coordinate system p,co,4>, centered on the contour center-see Figure 3.1. 
Note that it is not uncommon to use a left handed coordinate system. 




Figure 3.1. Coordinates. 



We can measure the components either by performing a Fourier analysis of the data from a set of 
individual coils measuring B n ((0), B T ((0), or we can construct integral coils which will do the job 

directly. For example, a "modified Rogowski coil", or "cosine coil", whose winding density 
(number of turns per unit length) nA = nocos(nro), each turn of area A, will give a signal which, 
when time integrated, is proportional only to X n : 
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d d\ I 

e=--\(B.n s dS) = --\\ (BJo))n A (co^Adco) } 



dt 



- - ^° An ° — 1 / [Tcos(n<») 1 + V /l n sin(n«) + // n cos(n<w)V|<i<W 

2^ d*[ J L v „ ' J] J 



2 J? 



3.3 



The elemental area dS = nAAdl, the unit length dl = aidco, and n s is the unit normal to the coil 
area. That is, the only contribution to the space integral comes from the term cos 2 (nco), because 
J q cos(ntt>)cos(mco)dco = % if m = n, otherwise = 0. If the winding density is proportional to 

sin(no)), the time integrated output is proportional to 8 n . To obtain the coefficients (i n and K n , we 
must wind a "saddle coil" with n w turns of width w varying as sin(no)) or cos(no)), so that for a 
"sin" saddle coil w((fl) = wq cos((fl), and 



e=--jU(B(co)n w w(a>)a l da>)\ 



dt 



3.4 



— _£JL 



M n ^ n w d 



dt 



In this case the elemental area dS = n w wdl = n w waid(0, the time integrated output provides the 
coefficient (X. Figure 3.2 shows a cosine coil which measures X\. Although it is not illustrated, a 

center return wound inside the Rogowski coil should be used. Figure 3.3 shows an unfolded "sin 
saddle coil" measuring Of course, we cold also use an array of coils placed on a contour, 
measuring independent B x and B n at different positions (different CO) and construct the required 
integrals. 



pitch changes sign 




'width' changes sign 




Figure 3.2. A modified Rogowski coil. 



Figure 3.3. A saddle coil. 
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Field components on a rectangle 



If we want to characterize the fields on a rectangular contour, we can make use of the fact that an 
arbitrary function in a plane can be expressed as 



3.5 



with c mp constant coefficients. Here we are working in a rectalinear coordinate system ^,T|, 
centered on the contour center, at R = R\, shown in Figure 3.4. 




Contour 1 



Figure 3.4. The geometry used in describing fields on a rectangle or square. 

On a one dimensional contour there will be degeneracy. Suppose we have a "modified 
Rogowski" coil whose winding density varies as some function fp(r|,£, ), so that the time 

integrated output is proportional to 

s p , T = jf p B,dl 3.6 

The subscript x refers to the tangential (normal) field component on the contour. We could also 
construct the signal s p x from individual measurements of B x around the contour. Further 

suppose that we express the tangential field itself in terms of our functions f as 

B T = T, c J m 3.7 
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Then we can write 

s =Yc {f f dl 3.8 

m i 

i.e. if we can calculate jTmfpdl for our chosen functions f, then we can express the coefficients c m 
through the measured parameters s P;T . In a similar way we can build a saddle coil of width f p 
whose time integrated output is then 

s =f/ Bdl 3.9 

I 

Expressing the normal field as 

B n = T,d m f m 3.10 

m 

we have 

s = V d if f dl 3.11 

p,n ZmU mV J mJ p 

m i 

Again, assuming we can calculate jjfmfpdl, the coefficients d m can be expressed in terms of the 
measured s p>n . 

We still have to choose the functions f p . One choice, which is used in 'multipole moments', 
discussed later, is pP, the p th power of a vector radius on the contour 1. This can be expressed in 
the form of a complex number as 

p=%+iri 3.12 

e.g. for p = 2 we have 

p 2 = e-V 2 + i{^) 3.13 

If the contour chosen is a square of half width and height a, then this form for the functions f 
gives 

ff m fndl = 4a {m+p) \ — - — I — — I if m and p are even 3.14 
, JmJp lm + l p + lj F 

j>f m f p dl = 0=0 otherwise 3.15 
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Then we would have the output from the 'modified Rogowski' coil 



s 



m 



-i for m, p even 



3.16 



= I> m 4a' 



m 



,(m+p) 



a a 

+ 

m+ 1 p + 1 



= otherwise 



3.17 



Equations 3.16 and 3.17 show specifically how, by including a finite number of terms (say p max 
= rn max = 5) we will end up with a set of linear equations relating the measured signals to the 
required constants c m . We must now solve them to obtain the coefficients c m as functions of the 
measured s PiX ; a similar procedure provides the d m as functions of the signals Sp 5n . The result is 
not as elegant as the Fourier analysis applicable on a circular contour, where a single coil can be 
wound to measure each individual Fourier coefficient, but I don't know another way to represent 
the fields on a square contour. Of course, instead of using these specially wound coils to 
measure s P;T and s P)I1 directly, the required integrals can always be constructed from individual 
coil signals of B T and B n around the contour 1. 

An example of a saddle coil for a particular f = T](l+^/Ri) is shown in Figure 3.5. Here R\ is the 
major radius of the contour center. These strange looking coils are actually useful for helping 
determine plasma position 




Figure 3.5. A saddle coil suitable for winding on a square vessel 
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4. PLASMA CURRENT 



Rogowski coil 



The plasma current is measured by a "Rogowski coil", which is a multi turn solenoid completely 
enclosing the current to be measured. Figure 4.1 shows an example, the placement of this coil 
around the plasma is shown in Figure 4.2. The transient plasma current generates a voltage £ 
which, for a uniform winding density of nA turns per unit length of area A, is (after applying 

Faraday's Law) 



from which I p is deduced after time integration. This is just a special case of our general model 
for how to measure the fields on a contour. Integration can be performed passively with a 
resistance-capacitance circuit, with active integrators, or numerically on a computer. In each case 
there is an associated 'integration time constant' Tint- The Rogowski coil must not be sensitive to 

other than the wanted field components, so that a center return must be used. The angle between 
the Rogowski coil and the enclosed current is irrelevant, as is the contour on which the coil is 
wound. 




4.1 





Figure 4.2. Coil 
placement around plasma 



Volts 



iRogowski 



Figure 4.1. A Rogowski coil 
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5. LOOP VOLTS, VOLTS per TURN, SURFACE VOLTAGE. 

Introduction 



The Loop Volts ei, also called the Volts per Turn or Surface Voltage, is used in calculating the 
Ohmic power input to the plasma. It also allows a calculation of the plasma resistivity Q p . £1 is 
also a useful measure of cleanliness: clean ohmic heated tokamaks usually have £1 ~ 1.5V. 

What we want to measure is the resistive voltage drop across or around a plasma. In a linear 
machine, this simply done by measuring the potential across the end electrodes with a resistive 
potential divider. A similar method can be used in a torus with a conducting vacuum vessel 
which has one or more insulating sections. In an all metal torus the voltage induced in a single 
turn pickup coil (a volts per turn loop) wound close to the plasma is used, as shown in Figure 5.1. 
However, the interpretation of the output signal is not trivial. Here I want to address two 
questions. The first is "What does a toroidal loop as shown in Figure 5.1 tell me?". The second 
is, "How do I measure the Ohmic power input into a plasma?". 

The single volts per turn loop 

The voltage across the toroidally wound volts per turn loop (subscript 1) is given by: 



Here subscript j refers to all fixed windings, such as the Ohmic heating, the vertical field, and the 
shaping winding. The plasma current contribution (subscript p) is brought out separately. We 
can arrange for the voltage of the loop to be measured with a high input resistance amplifier. 
Then Ii = 0, so that the first two terms on the RHS are zero, and 




5.1 




5.2 



If this signal is time integrated, then the result is exactly the poloidal flux x ¥, because 




5.3 
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Now consider the voltage £ p around the plasma. It is connected on itself (a torus) so that: 

e =0 = — (L 7 )+n 7 + Y — (M B ./.) 5.4 

p dt \ p,p ?) p p ^ dt \ P ,j j J 

Now remembering the definition of mutual inductance in terms of linked fluxes, we can always 
write the flux through circuit i due to current Ij in circuit j as the flux through another circuit k 
due to the current Ij in circuit j plus the incremental flux between the circuits k and i due to the 
current Ij in circuit j, A^k^j. Then 

M iJj = M k Jj + A<P M , = M.,7. + AV M } 5.5 
Then for example Mi i0h I oh = M P)0h Ioh + ^ v ,\-oh 
Thus we can write 

l df ^\ p,j jj df \ Up p) dt L^\ p,HjJ 

Substituting from Equation 5.4 gives 

£ ' = ~~dt ^ Lp,p Ip ^~ ^ pIp + ~dt ^ ' plasma -' oop ) ^ " 7 

where A^piasma.ioop is now the total flux between the loop and the plasma, provided by all 
circuits, including the plasma (plasma, ohmic heating, vertical field, shaping). If the plasma 
current is constant the volts per turn loop tells us the plasma resistance. A more elegant approach 
to seeing this is to use Poynting's theorem. 

Poynting's theorem 

Consider a number of non integrated flux loops, i.e. volts per turn loops, measuring d^/dt, all 
placed around the plasma on some contour 1, which might be the vacuum vessel. Figure 5.1 
shows the configuration. Note that the emf £ = -InKE^ will not necessarily be the same in each 

loop, because the contour 1 is not necessarily on a magnetic surface. 
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plasma 



Figure 5.1. A subset of volts per turn coils on a contour 1. 



To interpret what we are measuring, notice that the poloidal (subscript p) and toroidal fields are 
not coupled in Maxwell's equations, so that we can write Poynting's theorem for the poloidal 
fields alone. To remind you, the basic equations needed are 



VxB=//j 

VxE=-— 

dt 



5.8 
5.9 



Multiplying these by -E and B/(i respectively, adding, and writing the poloidal component, gives 
(<|) is the direction the long way around the plasma) 



d_ 

dt 



( D 2 \ 




f B ^ 




+ v« 


Ex^ 






v n ) 



5.10 



Integrating over the volume V defined by rotating the contour 1 in the § direction gives 
Lit 



dt 



+ 



$j^dV=$eB / ll 



5.11 
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Here we have used J V • (e x B p }iV = J(e, x B p )• dS^ = J 2nRE^B p dl . U is defined by (Ul p 2 )/2 



= j(Bp 2 /(2(io)dV. Note the integration is to the contour 1, not the plasma edge. Therefore the 
inductance is not just that internal to the plasma, which is usually called lj. Now use Ohms law 
j-B = C||E-B, and assuming IB^q-B^I « B^q gives = j^/oy. Therefore 



dt 



f L.r 

1 p 



■ 2 



+ \l±dV = I p (e) 

(J\\ 



5.12 



where 



(e) = -^-ieB T dl 



5.13 



What we find, from Equation 5.12, is that the ohmic input power \j§ 2 /o\\dY into the plasma must 
be evaluated knowing the poloidal distribution of both £ and B T around the contour 1, as well as 
the inductance Lj within that contour. For example, suppose the contour is a circle of radius ai, 
and 



B * = f^L 1+ ^ C0S( H 



5.14 



5.15 



Then we obtain 



5.16 



The inductance Lj is given approximately by (i.e. the "straight" circular tokamak ) 



r 



In 



+ 

\a p J 2 



5.17 



with a p , R p the plasma minor, major radius. The term (io/(47t)li is the inductance per unit length 

(toroidally) inside the plasma, and the In term represents the inductance between the plasma 
surface (r = a p ) and the contour 1 (r = ai). The approximation given is for a straight circular 

tokamak coaxial with a circular contour. That part of the inductance between the plasma and the 
contour is sensitive to plasma position. 
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Uses of the Volts per turn measurement 



We can deduce an average value of the plasma conductivity, <G>, by writing 



l^A = \i± dv=I{£) _l 



V 2 



5.18 



From this we can define a conductivity temperature T G . The conductivity deduced by Spitzer for 
Coulomb collisions is given by (there are corrections for the fact that, in a torus, trapped particles 
cannot carry current and so o must be reduced) 

Z^ln(A,) 



0-=1.9xlO 4 ^— — 5.19 



Then T G is defined as that temperature which gives a Spitzer conductivity (with Z e ff =1) equal to 
the average conductivity <o>, with an approximate value taken for ln(A s ). We can also derive an 
average "skin time", from the formula for the penetration of a field into a conductor of uniform 
conductivity <a>: 

^=- t ^ 5.20 



16 



The definition of energy confinement time for an ohmic heated plasma with major radius R p , 
cross sectional area Sphere we assume a circular minor radius a p so = 7ta p 2 ), total energy 
content W = 37tR p j pdS^, ohmic input power P h = I p 2 n p , can we written in terms of the 

poloidal beta value Pi = 87C/((loIp 2 )L pdS(K (discussed later) as 

E P nh 8Q p 16 
Combining equations 5.20 and 5.21 shows that 

— =A 5.22 

^ skin 

Therefore for ohmic heated plasmas, where typically (3i ~ 0.3, the currents penetrate 
approximately 3 times slower than the energy escapes from the plasma. 
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6. TOKAMAK EQUILIBRIA 

6.0. AN INTUITIVE DERIVATION OF TOKAMAK EQUILIBRIUM 
Introduction 

After having described how to measure the plasma current and loop voltage, the next most 
important parameter to measure is the plasma position. We will show how we determine both 
this, and certain integrals of the pressure and field across the plasma cross section (specifically (3i 
+ L/2), in section 7. The basic idea is that we want an expression for the fields outside the plasma 
in terms of plasma displacement A and (Pi + 1/2). We can only do this by knowing a solution to 
the plasma equilibrium, i.e. we must solve the Grad Shafranov equation. I deal with this later in 
this section, but first we can gain a physical picture of tokamak equilibrium by considering the 
various forces acting on a toroidal plasma. Also note that there are techniques to measure plasma 
position without recourse to equilibrium solutions, the so called "moments" method. However, 
the interpretation of this method (i.e. what has been measured) itself requires a knowledge of the 
equilibrium. 

The total energy of our system must be made up of 3 parts, 

W=W P + W? + W* 6.0.1 

where W p is the energy stored as pressure, Wi B is the energy stored in toroidal fields (poloidal 
currents), and W2 B is the energy stored in poloidal fields (toroidal currents). Once we have 

calculated expressions for these terms, we can obtain the required forces: the minor radial force 
F a = dW/da, and the major radial force Fr = BW/dR. By setting the net force = we will obtain 

the conditions necessary for equilibrium. We work with a circular cross sectioned plasma with 
major radius R and a minor radius a, and a/R « 1. Figure 6.0.1 shows the geometry. The 
poloidal coordinate is 0. We use <...> to mean an average over the volume (which is the same as 
an average over the plasma surface if a/R « 1) and V = 2jcR7ta 2 . A positive force is in the "R" 
or "a" direction (an expansion). 



45 



Magnetic fields and tokamak plasmas 



Alan Wootton 




The toroidal fields and pressures p 
present in a section of a toroidal 
plasma colum. the toroidal angle is 

4> 



Figure 6.0.1. Elemental volume discussed in deriving force balance. 

Energy associated with plasma pressure_w p 

The major radial force F P R exerted by the plasma pressure in expanding a distance dR is given by 
F P R dR = <p>dV, and the minor radial force F p a exerted by the plasma pressure in expanding a 
distance da is given by F p a da = <p>dV. The total energy WP is given by 

W = \pdV = {p)V = {p)l7t l Ra 2 : 6.0.2 

V 



so that 



r _ dW p _ 2V(p) 

da a 



_ dW p _ V(p) 
R dR R 



6.0.3 



6.0.4 



These forces were computed at constant pressure. 
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Energy associated with toroidal fields _Wi B 

The energy associated with poloidal currents is written as 

w B = hlL + hA + MII 6.0.5 

1 2 2 iii. 

Here Ii is the poloidal current in the plasma, and Ii e is the poloidal current in the toroidal field 
coil (subscript e for external). Ii is that poloidal current flowing in the plasma edge which 
produces a toroidal field equal to the difference between the internal toroidal field B^i and the 
external toroidal field B^e. By definition we have 

2 2ju 



Mo 



M \h l u = 6.0.7 



Now the circuits Ii and Ii e are perfectly coupled, so that Li = Mi The field Bi = (IoIi/(2tcR), and 

so 



M = A = "Ljf dV = Mq (r - (R 2 - a 2 / 2 ) - 6.0.8 
for skin currents. To get the forces we will need only the functional dependencies, namely 



_ 


iu 


da 


a 


dM l 




da 


a 






BR ~ 


R 


3M, 


—Mx 


dR ' 


R 



6.0.9 



6.0.10 



The forces will be computed at constant current. For example, the part of the force due to 
5/5R(LiIi 2 /2) is then written as (Ii 2 /2)e/eR(Li) = -(Ii 2 /2)L/R = -(Lili 2 /2)(1/R). Using Equation 
6.0.6 this becomes (Ii 2 /2)5/5R(Li) = <B ( | ) i-B() ) e) 2 >V/(2R(io)- Doing this for each component in 
Equation 6.0.5 gives 
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6.0.11 



6.0.12 



There will be no force from the term Li e Ii e 2 /2, because we keep the external currents constant. 

Energy associated with poloidal fields w? B 

Again we write the energy in the poloidal circuits as 

w!= h£ + LjL +Mil2he 6 .o.i3 

with I2 = I p the toroidal plasma current, and l2 e the poloidal currents in external windings. These 
windings are imagined to consist of a set which provides an external vertical field B z but induces 
no plasma current, and a set which drives the plasma current but produces no B z . We will need 

the radial derivatives (with respect to both major and minor radii) of the self and mutual 
inductances. All external currents will be kept fixed, so that 9(L2 e )/9a = 6(L2 e )/9R = 0. Only 

that winding producing a vertical field is imagined coupled through a mutual inductance to the 
plasma. 

We must make use of the inductance of a current loop: 




6.0.14 



The term \ accounts for the inductance between the center of the loop and the edge, at r = a. We 
also need the flux *P = M2I20 the external flux passing through the central aperture, which can be 
written as 

R 

^ = \l7tRB z dR 6.0.15 



Then we have 




6.0.16 
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dL, 


ju R 




da 


a 




dM 2 


d 




lltRB, 


dR ' 


dR 


v/ 2 ^ 


1 2e 



dM n 



da da \L 



2e 



— =0 



6.0.17 



6.0.18 



6.0.19 



Now we can deduce the forces: 



F s dW 2 



B 



da 



2V 
a 



( i?\ 



6.0.20 



2ft 



dw 2 B V 

dR ~ R 



Bi f, ( %R\ , n 2BB ff R 
—i In — -1+- + — 



6.0.21 



where 

B ea= ^k = tljL 6.0,22 

2 ^3 2 ^13 



Total forces 



We can now add Equations 6.0.3, 6.0.11 and 6.0.20 to get F a , and Equations 6.0.4, 6.0.12 and 
6.0.21 to getF R : 

a ^ 2// 2// 

F , . I [Af , _ , + iV + ^d^) +(p) j 6 . a24 

s /?!//„ V V a ; 2^ // a 2// j 



jU Q I 

with 5^ = L . Equilibrium requires F a = 0, so that 



(p) = | + v g L-gL/ + | 6 .0.25 
[ 2 A) 2 >"o i 
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Therefore if we can measure the difference in toroidal field with and without (when B^i = B^e) 
plasma present, we can measure the average plasma pressure. This is discussed later in dealing 
with "diamagnetism". 

The condition Fr = specifies the vertical field: 
aB a 



2R 



V a J 2 2 B 2 



da 



6.0.26 



This was the field we used in section 1, Field lines and flux surfaces, to plot out the flux surfaces 
which result from a combined circular filament and a vertical field. 

6.1. THE FLUX OUTSIDE A CIRCULAR TOKAMAK 

Later we will use the expression for the flux outside a circular tokamak. It can be considered to 
come from two sources, that from the external maintaining fields \|/ ex t and that from the plasma 
itself, V|/ p . In the previous section we derived an expression for the vertical field B z necessary to 

maintain a circular equilibrium (Equation 6.0.26). While the major radial term appearing as 
f 8i?^ 

lni — in Equation 6.0.26 clearly refers to the geometric center R„ (it comes from the 

^ a y 

inductance of a plasma with radius a with a geometric center R g ), it is not obvious to what radius 
the term outside the square brackets refers. It could be either R g , or the coordinate itself, so that 
B z oc 1/R. In the former case, in a right-handed cylindrical coordinate system (R,(|),z), the flux 
would be derived from y R 2 , or in our local coordinate system (p,co,cf)) based on R g (See 
Figure 1.7) 

y/ = £ + fcj cos(ftj) + & 2 cos(2<y) 6.1.1 

with ki a constant. The constant is unimportant, but the cos(2a)) term means that such an external 

field would introduce ellipticity, and we have specifically considered a circular plasma. 
Therefore we must take B z oc 1/R, and in the coordinate system (R,(|),z) this is derived from a flux 

uJ n R\ (%R\ 1 
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where A = /?, + ^ - 1 . In the local coordinate system (p, G), (])) the necessary B z is then derived 
from a flux (ignoring the constant of integration) 



Y, 



<u Q I p pcos(G))[^ (SR\ 



4/r 



+ A-0.5 



1 



6.1.3 




current filament center R 



geometric center R g 



Figure 6.1.1. The geometry used in relating the geometric and current filament centers 

Next we come to the flux \|/ p produced by the plasma. Outside the plasma, where there is no 
current and no pressure, the fields and fluxes we are looking for must be able to be constructed 
from those due to circular filaments. This will not be true inside the plasma. For a first 
approximation we will model the flux \|/ p as being due to a single circular filament with current 
Ip. If we position this filament at a position R c then in a coordinate system (p c , (0 C , (])) based on 
the filament we have shown that the flux is well represented by 



Y„~ 



In 



In 



-2 



Pc. 



CO; 



'4^ 



-1 



1 



6.1.4 



However, we have to decide where to place this filament, that is, what to choose for R c , and how 
p, Pc, 0) and G) c are related. Because we have derived \|/ ext for a circular cross section (in the 
calculation of B z we used inductances for circular current path), we must place the filament so 
that, in the coordinate system (p, (0, (])) the surface \|/total = Vext + \|/ p = constant is a circle of 
radius at p = a. 
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Suppose we place the filament at a position Ro = R c + A. From Figure 6.1.1, and expanding in 
the small parameter A/r, we then derive that 



p) =/7 2 +A 2 -2A/7cos(180-<w) 
p 2 = p] + A 2 - 2 Ap c cos(co c ) 



6.1.5 
6.1.6 



from which we have 



Pc~P 



1 +— cos(<») 

p 



6.1.7 



and 



cos(<» c ) ~ cos(co)(l - — cos(co)j + — 



6.1.8 



Substituting for p c and G) c in terms of p and (0, and using Ro = R c + A, we obtain an expression 
for the total flux V|/ to tal (P> <Q, (])) in the coordinate system based on the geometric center. Keeping 
terms of order A/r and cos((0) only (i.e. neglecting elliptic distortions to any surface) we find 



Y, 



total 



In L v p J J An L ^a J 2 p 2 J 



6.1.9 



To ensure a circular outer contour (at r = a) we set the cos((fl) term to zero, that is we set 



a 2R g 



a ^Ax 1 

i A + - 

V 7 



6.1.10 



Substituting for A form Equation 6.1.10 into Equation 6.1.9 gives the final expression for the flux 
outside a circular tokamak: 



Y, 



total 



271 



1 



-2i + 



An 



ln laJH A+ 2 



1- 



a ' 1 



\pj ) 



6.1.11 
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6. CIRCULAR EQUILIBRIUM 

Derivation of the Grad Shafranov equation 

Most analysis of magnetic measurements rely on an understanding of equilibrium. The most 
general cases are derived from integral equations, which do not constrain the plasma shape, and 
are discussed later. However, an insight can often be gained by considering the case of a circular, 
low beta, large aspect ratio (a/R p «l) system. For this we must solve the so called Grad- 

Shafranov equation. The basic equation for equilibrium is that of pressure balance: 

jxB = V/? 6.1 

from which B Vp = 0, i.e. there is no pressure gradient along a field line; magnetic flux surfaces 
are surfaces of constant pressure. Also from 6.1 we have j-Vp = 0, so that current lines lie in a 
magnetic surface. Equation 6.1, together with Maxwell's equations, is all we need. 

The conditions for the applicability of equation 6.1 are interesting. It is NOT correct that it only 
applies to the ideal MHD plasma, which would imply that it is incompatible with diffusive 
phenomena. In fact it is sufficient that 

1) The plasma directed kinetic energy is much smaller than the thermal energy. Inertia is 
ignored. 

2) The gas kinetic pressure is nearly isotropic. This is still possible for a collisionless 
plasma (k » 2tcR) as long as the times considered are long (t » X/v t h, i.e. t » x ee , T e i with x ee , 
x e i the electron -electron and ion-ion collision times). 

3) The local distribution function deviates slightly from a Maxwellian even for a straight 
system because of the finite Larmor radius ru 8f ~ rL-5f(r)/5r. This leads to a deviation from 
pressure isotropy by an amount 8p ~ pri/a, with a the minor radius. In toroidal geometry the 
particles are displaced from a magnetic surfaced by amounts d ~ rLq for passing particles and d ~ 
r Lq£~ 1/2 f° r trapped particles (£ is the aspect ratio r/R, q is the safety factor). Then there is a 
deviation from a local Maxwellian distribution function f by an amount 8f/f ~ d/a, and an 
associated pressure anisotropy 8p/p ~ rLq/a. Therefore we should have rLq/a « 1. i.e. the 

plasma can be considered as a gas of Larmor circles produced by charges rotating in a magnetic 
field. 
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Figure 6.1. Nested flux surfaces. 



Figure 6.2. Field lines in a surface. 



We introduce \|/, the poloidal flux per radian in (]). This is proportional to the poloidal flux within 
each surface; it is constant on a flux surface. Figure 6.1 shows an example of nested flux 
surfaces, and Figure 6.2 shows field lines lying in a flux surface (followed 10 times in the 
toroidal direction). We will show that the current lies in a flux surface, but current lines do not 

follow field lines. In Figure 6.1 the each surface is described by a radius r = r + ^ A„ cos(na>), 

where we have also shifted the origin of each surface with respect to each other. The example 
shown has a circular innermost surface, a mostly elliptic (n = 2) central surface, and an outer 
surface with a combination of ellipticity (n = 2) and triangularity (n = 3). Note that \j/ = RA^, 
where A is the vector potential (note we sometimes use the total flux m = 2n\\f). In the 
cylindrical coordinates (R,(]),z) of Figure 1.7 the magnetic fields are derived from 



B R = 



R dz 



6.2 



z R dR 



6.3 



which satisfy VB = 0. We also introduce a current flux function f, which is related to the 
poloidal current density so that 



R Rdz 



6.4 
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h = — 6.5 

Jz R dR 

Now Amperes law (ifj = V x B gives 
Mo dz 

• 1 B ( RB ^ 6.7 



z jU R dR 

so that, comparing Equation 6.4 with Equation 6.6, and Equation 6.5 with Equation 6.7, we have 
fi f=RB, 6.8 

That is, the function f includes the total current in the windings producing the toroidal field. 
Since from Equation 6.1. j-Vp = 0, using Equations 6.4 and 6.5 for j gives (9f/6R)(9p/3z) - 
(d&dz)(dp/dR) = 0, or 

VfxVp = 6.9 
Since p is a function of \\f, i.e. p = p(\|/), we must have f = f(\|/) as well. 
Now we want to derive the basic equilibrium equation in terms of \|/. Write Equation 6. 1 as 

where subscript p means poloidal and is a unit vector in the § direction. Now Equations 6.2. 
and 6.3 can be written as 

*,=\frr**>) 6.11 
and Equations 6.4 and 6.5 can be written as 

j,4( v/xe ^) 6 - 12 

Substituting 6.1 1. and 6.12 into 6.10, (remember e^-Vy = e^-Vf = 0) gives 
B i 

l Vf + ^Vy/=Vp 6.13 

R R 
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Writing Vf(\|/) = (df/d\|/)V\|/ and Vp(\|/) = (dp/d\|/)V\|/, and substituting into 6.13, gives 

6.14 



„ dp „ df 
= R—+B, 



dy/ v dy/ 



or 



= R d P , y d f 

dy/ R dy/ 



6.15 



The equilibrium conditions restrict the possible current distributions in the plasma ; instead of 
a two dimensional distribution j(])(R,z) a one dimensional one is obtained, depending on the two 
profiles dp/d\|/ and fdf/d\|/. 

We can write in terms of \|/ using Amperes law, (igj = V x B. Substituting from Equations 6.2 
and 6.3 for the components of B in the toroidal component gives 



6.16 



Substituting for from Equation 6.16 into Equation 6.15 finally gives the Grad Shafranov 
equation: 



R— 



d fl dy/) a> 



dp 



df 



dR^RdRJ + dz 2 " M ° R dy/ Mof dy/ 



6.17 



In the local coordinates (0,r,<f>) of Figure 6.3, based on the plasma major radius (the geometric 
center of the outermost circular surface) R g this equation becomes: 



1 3 ( d\ | 1 3 2 
r 3rv dr) r 2 d0 2 



¥- 



1 



(R g +rcos(0)) 



cos(0) 



a sin(0) a 

dr r dO. 



¥ 



6.18 



Note the ordering (0,r,<f>) in which the system is right handed. The poloidal angle is measured 
from the outer equator, so that x = (R + rcos(#))cos(^), y = {R+ r cos(6>))sin(^), and 
z = rsin(#). The metric coefficients are h e = r , h r = 1 and h^ = R + r cos(#) . If a new variable 
01 = 71 - is defined, so 0i is measured from the inner equator, then the system (r,9i,<f)) is also 
right handed. We have mostly used this latter system. 
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Figure 6.3. Geometry of the local coordinate system used. 



Solving the Grad Shafranov equation 

In a large aspect ratio expansion let us write 
y/= Wo(r)+ Vi(r,0) 



6.19 



where the first term describes a set of concentric circular surfaces, and the second term is the first 
order correction. Substituting into the general equation (Equation 6.18) gives the zeroeth or 
leading order equation for 



r dr^ dr 



dy/ 



dy/ 



6.20 



The first order equation for \\f \ is rearranged as: 

l~l d ( d\ 1 d 2 1 cos(e)dy 
-— i r— +-t — r ¥1 — 



R„ dr 



d_ 

dr 



— — - VI/, 

d\\f J d\|/ 



-2(i o R g rcos(0) 



dp(¥o) 
d\|/ 



6.21 



Suppose each surface \|/ is displaced a distance Ai(\|/n(r)), so we can write \|/ as 
¥= Wo + Wi = Wo -\( r )^r = Wo-\(r)cos(0)^- 



dR 



dr 



6.22 
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Substituting this form for \\f\ (= -Ai(r)cos(0)d\|/o/c3r) into the first order equation, Equation 6.21, 
gives: 

dr^ Be0 drJ-R g ^^ r dr 623 
1 dw n 

where B an = has been used. 

R g dr 

The solution of this equation gives the horizontal shift Ai of a flux surface away from the 
geometric axis, as illustrated in Figure 6.4. This, together with the solution of Equation. 6.22, 
gives us what we want. If we wanted to consider non circular surfaces, we could write 

r = r o + ^A n (r)cos(n0) for each surface. Then \|/(r) is a constant (by definition), so that 




Figure 6.4. Displaced circular flux surfaces. 
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The poloidal field at the plasma edge 



To continue, we need an expression for the boundary field Be(a) from the plasma equilibrium. 
We will ultimately use this to match the external and internal solutions at the plasma boundary r 
= a, where B T (= Be for a circular plasma) must be continuous. Now 



1 



dif/ 



1 dif/ 

R dr ~ (R g + r cos6>) dr 



6.24 



Using Equation 6.22, with Ai(a) = 0, gives 



BJa) = B 0O (a) 



r 



i- 



a ( dA 1 



— + 



dr 



cos(#) 



a J 



6.25 



with Beo(a) the zero order poloidal field at the boundary. Now we need to find an expression for 
(dAi/dr) a in Equation 6.25. Integrating Equation 6.23 (with the condition dAi/dr = at r = 0) 

gives the displacement of the magnetic surfaces for the zero order pressure distribution po(r) and 
zero order poloidal field Be (r): 



JA, _ 2jU Q 
dr ~ rR.lt 



2 

r Po 



00 V 



-h Po 



+■ 



B: 



00 



2ju , 



\rdr 



6.26 



By defining (we will be more careful with definitions later) 



Mpsi aB an (a)i 



a 2 B 2 S0 (a) ' 



6.27 



for our circular case, and 



/.■=■ 



Mo 



— f 



Bt 



^dV=- 



Bl 



2ju " aB 2 eo (a)\ \2ju 



-rdr 



6.28 



where B p is the poloidal field, Sp^ the plasma cross section in a poloidal plane, and V p the 
plasma volume, for our circular case we obtain 



fdAi 
^ dr 



6.29 



Now we can substitute Equation 6.29 into 6.25 and obtain 
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B d (a) = B so (a) 1 + — A cos(0)^ 



R 



6.30 



where 



A=J3 I + 1 - L -1 



6.31 



That is, we now know the distribution of poloidal field at the boundary of a circular plasma 
including terms of order (a/R). 

Simple current distributions 

Consider a case where the zero order (circular cylinder) current density is given by 



J(f> ~ J(f>0 



'>o(f-* 2 ) 



6.32 



with < x = r/a < 1 representing the normalized minor radius, 
l 

0.8 
0.6 
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Figure 6.5 The normalized current density distributions j^/j^O). 

Figure 6.5 shows these profiles for integer a = through 5 as a function of normalized radius x. 
The total current J is given by 



Ka 
1+ a 



J</>0 



6.33 



so that we can write a normalized current density j r as 
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j r = —f L =(l + a){l-x z ) 



6.34 



Figure 6.6 shows j r as a function of x for a between and 5; these curves represent the shape of 
the current density profiles with the constraint that the total current J is constant. 




0.2 0.4 0.6 o.s l 

x 

Figure 6.6 The normalized current density distributions j r (x) 



The poloidal field Be is given by 

M-**r) 



B e ~ BeoO-)' 



6.35 



x 



with Bqo(1) the value at r = a, or x = 1: Bqo(1) = |l{)J/(27ia). The ratio Be/Beo(l) is shown in 
Figure 6.7 as a function of x for a between and 5. 

5_ 




Figure 6.7. The normalized poloidal field distributions Be/Beo(l) 
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The "self inductance per unit length" is expressed as 



6.36 



The solution is written in terms of PolyGamma and EulerGamma functions: 
lj = EulerGamma + 2 PolyGamma[0, 2 + a] - PolyGamma[0, 3 + 2 a] 

However, a polynomial fit in a can also be used: 

h = 0.509619 + 0.462798a- 0.0630876a 2 + 0. 00443746a 3 



6.37 



6.38 



Figure 6.8 shows the value of lj as a function of a for both the exact solution (Equation 6.37) 
and the polynomial fit (Equation 6.38). The fact that you cannot distinguish the two lines 
demonstrates that the fit is good. 




a 

Figure 6.8. The value of lj as a function of a. 

We must determine how to choose the free parameter a. Two ways are suggested. First, if the 
safety factor (discussed later, but q = (rB^RBe)) q(l) at x = 1, (i.e. at r = a) and q(0) at x = 0, 
(i.e. at r = 0) are known, then 



q(0) 



6.39 
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We must choose q(0). Figure 6 .9 shows the resulting a as a function of q(l) for two assumed 
values of q(0): the long broken line is for q(0) = 0.8, and the short broken line is for q(0) = 1.0. 
An alternate prescription is to choose a so that the position xi of the q = 1 surface is in 
approximately the correct position. For example, suppose that xi = l/q(l). Then the equation for 
local safety factor 



q(x) = q(\) 



x 



(i-(i-*r) 



6.40 



can be solved to give 



a = 



9(1)' 



6.41 



Figure 6 .9 (solid line) shows the resulting a as a function of q(l). 




Figure 6.9. The parameter a as a function of q(l), chosen such that 1) solid line: 
xi = l/q(l), 2) long dashed line q(0) = 0.8, and 3) short dashed line q(0) = 1.0. 



We have now uniquely determined a in terms of q(l), and as such there is a unique value of q(0) 
for each q(a). This is shown in Figure 6.10 
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Figure 6.10. The value of q(0) as a function of q(l), with the constraint that xi = l/q(a). 

For completeness, Figure 6.11 shows the various normalized q(x)/q(l) profiles for various values 
of a. 



0.8 
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Figure 6.11. Normalized profiles of q for various values of a. 

The surface displacements: the Shafranov Shift 

If we want to obtain actual values of Ai(r), we have to assume functional forms for the current 
and pressure distributions, and integrate equation 6.26. We have already discussed a simple 
current distribution; we must add a simple pressure distribution. Let the zero order pressure be 
given as 

pW = p (0)(l-/) r 



6.32 



so that the poloidal beta value, by definition, is 
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Pi 



2^oPq(0) 

(1+ rKo(i) 



6.42 



Now we can use the normalization x = r/a in equation 6.26 to obtain an expression for the surface 
displacement: 



2x 



W+iy(i--7-A(i-(i- 2 r>Ki-(i-- 2 rj-i 

x J 



6.43 



Where £ = a/R g . The Shafranov Shift is defined as the distance between the magnetic and 
geometric axis: A s = Ai(l). Using a power series expansion for Ai(x) up to terms x 6 a general 
expression can be derived: 

1 



2a a 2 
I 9 72 



A =— I 4 



+ 



6.44 



\ PMy±}) \ 2 , iig , 5a 2 | (r-i)(r-6) Mh) | I 

[ 4(1 + a) 2 1 9 18 6 9 J J 



For the simple case of a flat current profile (a = 0, lj = 0.5) and a parabolic pressure profile (y : 
1) we obtain 



6.45 



Matching vacuum and plasma solutions 



Returning to the field outside the plasma, we must match the vacuum field to the solution for 
Be(a). The vacuum field is given by the solution to (VxB)^ = 0. Expressing Br and B z in terms 
of \|/ (Equations 6.2 and 6.3), (VxB)^ has the form of the LHS of Equation 6.18. This has a 
solution of the form (for r « R g ); remember we are outside the plasma: 



f A 



¥ = 27cy(r,e) = ^ R g I p 



V V r J' 2 



In 



+ ■ 



In 



8R„^ ^ 
r J 



Ci ^ . 
+ — + C 2 r cos(0) 
r 



6.46 
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where the constants Ci and C2 are given by the boundary conditions (B-n) = at the plasma 
surface (B r (a) = 0) and the continuity of tangential field at the plasma surface. At infinity \\f takes 
the form C2rcos(0), i.e. this represents the part of \|/ from the external sources. The first term on 
the RHS is the flux from a filament. Substituting 6.46 into 6.24 gives an expression for Be(a), 
which is matched to the plasma solution equation 6.30 by taking 



"T - C 2 = lni — L + 2 A 

a " v a J 



6.47 



-1 dw 

Now since B r = — the requirement B r (a) = means that the coefficient of cos(0) in 6.46 is 

rR g dO 

zero at r = a. Therefore we obtain 



q = a 2 (A + 0.5) 

C 2 = -ln^-^j + A -0.5 



6.48 
6.49 



Thus we finally obtain an expression for the flux outside the plasma r > a: 



2k 



XTJ'V 



6.50 



The field components are given by 
Be -R3r' Br -"R CT r 30 



6.51 



i.e. 



rr 



2\ 



Be(r ,e) = -^-^ i4 



2%r 2ttR ( 



V 



J 



(A + i| + ln 



r ^ 1 

r 



V d p/ 



cos(0) 



6.52 



B r (r,e>=-^ 



4tcR 



1 --2- 

V r J 



(A + i| + ln 



( \ 



sin(e) 



6.53 
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We can also separate the total flux \\f into external (\|/ e ) and plasma (\|/ p ) parts. At large r, the 
expression (ln(8R g /r) -2) really approximates an expression which goes to zero. Therefore at 
large r we have 



ln^ — L j +A-—jr cos(6>) 6.54 



leaving 



V 



MortJ , (MA , faVf A nl 6 ' 55 
— — ^ — ^-i —In, — + i A + — | icos(^) 



2n 

i.e. the external field at the geometric center (R = R g ) has a component 



5 =• 



lni — L +A--, 6.56 

V a J 2) 



This is the component of the field which must be produced by external conductors to maintain a 
tokamak equilibrium. 

Figure 6.12 shows contours of the poloidal flux for a plasma with geometric center 0.06 m 
outside the coordinate center, minor radius a p = 0.265 m, and A = 2. At the inner equator an X 

point appears, where d\|//9R = 0. At this point the external vertical field and the poloidal field 
from the plasma cancel, so that B z (= -Be) = 0. The flux surface through this point is called the 
separatrix. This X point approaches the plasma surface as (3i increases. As this happens the 

surfaces are strongly distorted from circular, and the above equations starts break down. Figure 
6.13 shows the poloidal flux \\f = v J / /(2tc) in the plane z = for a plasma with minor radius a p = 
0.2m, a major geometric radius R g = lm, current I = l/\io, and A = 2. The external source 

(vertical field) and plasma components are shown separately. Figure 6.14 shows the minor radial 
(B r ) and poloidal (Be) field components for this plasma, on a circular contour concentric with the 

plasma minor radius. 
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Figure 6.12. Flux contours for a (badly drawn) circular plasma with geometric center 
0.06 m outside the coordinate center, minor radius 0.265m, A = 2. 
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Figure 6.13. The external poloidal flux \|/ = v J / /(2tc) in the plane z = for a plasma with minor 
radius a p = 0.2m, a major geometric radius R g = lm, current I = l/(io, and A = 2. 
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Figure 6.14. The field components B r and Be on a contour with minor radius 0.3m 
placed outside and concentric with the plasma described in Figure 6.13. 



More complicated configurations. 



It is easy to make non circular cross sectioned plasmas. Adding a quadrupole field produces an 
elliptic deformation. Adding a hexapole field produces a triangular deformation. There are some 
analytic expressions available for non circular cross sections, but in general numerical solutions 
to the basic Grad Shafranov equation must be used. 
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7. Position and (3i + \J2 for the circular equilibrium 

An 'exact' circular equilibrium 

Consider the pickup coils shown in Figure 7.1. Two coils (B^i, B ro 2) measure the poloidal field 
at the inner and outer equator, each at a distance b from some center at R = Ri. A saddle coil 
measures the difference in flux between the inner and outer equator. Instead of the partial flux 
loops (i.e. the saddle coil) shown in Figure. 7.1, two complete flux loops could be used, as in the 
plan view of Figure 7.2. We want to know what these coils tell us for our circular equilibrium. 




Fig 7.1. Saddle and Bg) coils Fig 7.2. A plan view of poloidal flux loops 



The original equilibrium was derived in a coordinate system (r,9,<f>) based on the plasma 
geometric center R g . We must translate the results to a coordinate system (p,CO,(|)) based on for 
example the center of a circular contour centered at R = Ri, which might be the vacuum vessel 
center. Details are seen in Figure 7.3. 
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sma with geometric center R 

- R 



g 



contour with center R 



Figure 7.3. The coordinate systems associated with the geometric center R g , 
and an external contour center R\. 

Then R = R\ + pcos((fl), z = psin((fl), and let A g = R g - R\ be the displacement of the plasma 

geometric center from the contour center. Remembering that we are only dealing with r > a, the 
plasma minor radius, and taking Ag/a p « 1, we can rewrite Equation 6.50 as 



In In 



. v p J ) 



An 



2\ 



1-H 



A+^-l+ln 



f \ 
P 



2#,A„ 



cos(<») 
7.1 



The tangential and normal field components on the circular contour are then 



a ^ 2np 4nR l 

M I P ( a- 



f 2\ 

1 + - JL 

V. P J 



A + - | -f In 



P_ 

\ a P j 



-l+- 



27?, A „ 



BAp,0)) = - 



4nR, 



2\ 



I — 

V P'J 



A +- | + In 



f \ 
P 



P 
7.2 

sin(<») 



cos(<») 



7.3 



Now turning to the coil configurations of Figure 7.1 and 7.2, the outer coils at p = b, G) = 
measure B^, and the inner coils at p = b, u) = n measure B ro i. The saddle or flux loops measure 
the difference in flux between the inner (^i) and outer (¥2) equator. Define 

then 
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and 



- B. 



oil 



art 



+ B ± = 



In 



\ a P j 



7.5 



b 



2R t 



f 



In 



+ 0.5 



f 2^ 

a 

1 e. 



V 



+ • 



^7 



2 



^ 2 a 



V 



2\ 



1+- 



7.6 



Therefore, assuming circular equilibrium, we can obtain the geometric displacement of the 
plasma from the vessel center, and the parameter A = (3r + L/2-1 if we have an estimate of the 
plasma radius a p . To get a p we must iterate our Equations 7.5 and 7.6, starting with an assumed 
value (say the limiter radius aii m ), and replacing it at step n by the maximum possible radius 
defined by A g and the specified limiter geometry at step (n-1). Such a procedure is shown below 



Iterations to determine magnetic 
configuration from fields. 



1. Set initial guess for plasma minor 
radius ap = limiter radius. Measure 
poloidal fields. 



2. CalculateA andA using 
estimated ap. 



3. Calculate plasma geometric center from 
estimated ap,A and A 



4. Estimate new radius ap 
from plasma geometric center 
and limiter locations. 



A more sophisticated method is to use a modified Rogowski coil (nA cos((0)) to measure the 
part of B^p,©) proportional to cos((fl), and a saddle coil (width oc sin(G))) to measure that part of 
Bp(p,Cfl) proportional to sin(G)). Such coils were shown in Figures 3.2 and 3.3. The coefficients 
measured (after time integration) will be, for a coil on a circular contour of radius p, 



p 



2R, 



f 2\ 

1 + — L 

V P J 



A + ^| + ln 



f \ 
P 



1+- 



27?, A „ 



7.7 
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Mi = 



P 

2R, 



2\ 



P 2 ) 



A + -J +ln 



f \ 
P_ 

\ a P j 



7.8 



From these two equations we obtain 



£(4-A)+-£r = A, + ^-(A+0.5) 



4R, 



2R, 



7.9 



— +//!>= A + log 



A, 



7.10 



from which we can obtain the two unknowns A and A g , assuming a p (which must be determined, 
together with A and A g , by iteration). Note that the term ln(r/a) represents that part of the 
inductance between the plasma surface and the coils. 



Extension of position measurement to non circular shapes 



The previous section showed how, for our circular equilibrium model, the geometric center could 
be measured. We would like to generalize the method to allow non circular shapes. In fact, the 
exact same diagnostics can do this, with certain important assumptions about the major radial 
dependence of the flux x ¥. 




Figure 7.4. Flux loops and poloidal field pick-up coils outside a plasma. There is a 
coordinate system based on the center R = R\ of the two measuring coils and flux loops 
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(sensors). The sensors are at a distance b from R. The plasma geometric center is 
displaced A g from Ri. 

Consider the situation in Figure 7.4. The measurements of flux and field are supposed to be at 
the same location, either at point 1 (shown at the smaller major radius) or point 2 (larger major 
radius). If the loops have the same value of flux, then they sit on a flux surface. Using the saddle 
coil as shown in Figure 7.1, if the saddle coil output is zero then the inner and outer toroidal legs 
are on a flux surface. We expand the flux as 

vp = vi/ + a 7 11 

plasma loop ()R 

with A the distance between the loop and the plasma surface. This expansion goes badly wrong if 
there is an X point between the measuring loops and coils and the plasma surface, because in this 
case there is a point between the measuring loops and the plasma where d^/dR = 0. Since = 

l/(27iR)9 v I / /9p, we can write 

¥1=^ + 2*^ A, 7.12 

* 2 =^+2^AA 2 7.13 

Now A2 = b-a p -Ag, Ai = b-a p +A g , Ri = R\-b, R2 = Ri+b, so that Equations 7.12 and 7.13 give 



2^(11, -b) + B m2 (R l+ b)] 



-(b-a p ) 7.14 



Therefore we can measure A g if we know the plasma minor radius a p . In limiter geometry a p is 
determined by A g and the assumed known limiter geometry, so that as described before iterations 
are necessary to provide both a p and A g . If necessary terms in 5 2v I7dR 2 can be measured using 
more B a coils (two at the inner, and two at the outer, equator). If the plasma is bounded by a 
separatrix, then the whole expansion is useless anyway, and we must turn to integral relationships 
discussed later. 

Extension of fa + 1/2 measurement to non circular shapes 

From our simple circular equilibrium described in section 6, we know that the external field 
required to maintain a circular low beta plasma is given by Equation 6.49: 



B Z =- 



In 



P 



+ A-± 

2 



7.15 
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Therefore if we know all the currents in the external conductors, the plasma current, minor radius 
and major radius, we can calculate Bz at the geometric axis R = R g , and obtain a value for A = (3i 
+ lj/2 -1. The result (Equation 7.15) has been found by numerical simulation to be accurate even 

for non circular discharges. Therefore if we can calculate the maintaining field, and know the 
plasma minor radius and geometric center, we can calculate Pi + lj/2. It becomes difficult to 

calculate the maintaining field if we have an iron core present; we will discuss this later. 
However, there are more general techniques which do not explicitly require this calculation of 
fields from currents, and we also discuss these later. 

Non-circular contours. 

So far we have only considered poloidal field measurements which can be made on a circular 
contour. In many cases the simplest contour to use is that of the vacuum vessel, which is often 
non circular. What do we do then? In section 3 we discussed how fields on a rectangle can be 
characterized, and suggested that coils wound to measure "moments" might be useful. 
Specifically, we discussed modified Rogowski coils which would measure 

s Ptt = $f,BjU 7.16 

i 

and saddle coils which would measure 

s =LfBdl 7.17 

p,n J J p n 
I 

where the functions f p were taken to be the p th power of a vector radius on the contour 1. The 
geometry is re drawn in Figure 7.5. 
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Figure 7.5. Geometry used for a non- circular contour. 



An analogy with our cosinusoidally wound Rogowski coil and sinusoidally wound saddle coil 
would be a modified Rogowski coil measuring the first symmetric (in vertical position) moment 



7.18 



and a saddle coil measuring the first asymmetric (in vertical position) moment 



p 



2^ 



1 + — £ 

V P ) 



A + -| + ln 



( \ 
P_ 

\ a P j 



2ft A 

1+^T^ 



7.19 



To interpret what these coils will measure, we can write an equation for the components B T and 
B n on our chosen contour (for a rectangular contour they will be either or B^). Because the 

only variables in the equations for an assumed circular equilibrium are the geometric 
displacement A g , A and minor radius a p , we must be able to derive expressions for the measured 
coil outputs si )T and si >n in terms of these variables. For example, if our contour is a square of 
half height and half width a, centered at R = Ri, we must find 



S l,r = ' S l,X / P' A ' A g' a P ' a '^) 



7.20 



7.21 
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In principle we can then solve Equations 7.20 and 7.21 to give expressions for the required A 
and A g . Iterations will be necessary because we will find a p entering the final results, which itself 
is only determined once A g is known. 

It may be that the particular geometry of the contour precludes performing the analytic integrals 
(Equations 7.18 and 7.19, using the analytical representation for the equilibrium fields outside a 
circular plasma). In this case the problem must be solved numerically. Equations 7.18 and 7.19 
are solved numerically for many different A g , A, and a p (we assume the contour geometry does 

not change!) A regression analysis is then performed, in which solutions of the form 




7.22 




a .contour geometry 



1.23 



J 



are sought. 
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8. SOME FUNDAMENTAL RELATIONS 

Geometry 

In section 6 we derived an analytic expression for the flux outside a large aspect ratio (a/R g « 
1) circular plasma, which we used in section 7 to interpret magnetic field measurements. Here 
we want to derive some relationships which are of general use: we shall test them for the analytic 
equilibrium of section 6. Use a right handed cylindrical coordinate system (R,(|),z), usually 
symmetric w.r.t. rotations around R = 0. Referring to Figure 8.1, there is an axisymmetric region 
V, which completely encloses the plasma. The cross section of V in the poloidal half plane (<|> = 
0, R > 0) is S§, and the boundaries are S n and 1. dV is the volume element on V, dS^ the area 
element on S^, dS n is the area element on S n and dl the line element on 1. Therefore 

dV = 27lRdS^ 8.1 

dS n = 27tRdl 8.2 



z 




Figure 8.1. Geometry 



Normal and tangential derivatives on 1 are d/dn and d/dx. The positive orientation on 1 is that 

lies on the RHS. Note this is DIFFERENT from what was assumed in section 6. The general 
question is: "How can we derive information on the magnetic field and plasma in V (or S^) from 

the fields measured on 1?" 
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Field representation 



We shall try and be as general as possible. First we quickly re-derive the equilibrium equation. 
For the case of axisymmetry, we discussed in section 6 the description of the fields using the 
form 



2k 2k 



\ v p J 



4k 
8.3 



f 2\ 

V p 2 j 



A+|| + ln 



f \ 
P 



cos(<») 



where 



F=RB,=pf 



8.4 



and \|/= RA(j) 



From VxH = j, we can also write 



1 



j=--L yN(/) + V - \xV</> 
p \pJ 



8.5 



where the operator L* is defined for axisymmetric scalar fields as 



Ly/ = pFV 



pR 2 ) 



8.6 



The operator L* satisfies the equation 

i\j/=-pRj $ 

For uniform permeability (J, = (lo we have L* reducing to the operator A* 

2 ( VmA d ( 1 dyA d 2 w 
A y/= R V •[ — j- = R— +—7- A 

Y V R 2 J dRKRdRj dz 



8.7. 



8.8 



That is, A* is the operator on the LHS of the Grad Shafranov equation. We also introduce the 
operator L: 



Ly/=-V»{pVy/) 

P 



8.9 
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which reduces to the Laplacian A for uniform permeability. In a current free region we can use 
the representation H = Vg, and V B = is then equivalent to Lg = 0. 

Identities 

Now we turn to some identities. Green's first identity for L* is: 

\\yrt%dS $ = l\yf^dl- \—VyV®dS, 8.10 

Green's second identity (Green's theorem) is: 

f J-(^0-0L>W =L± f w ^.- @ ^) d i g.ii 

Both of these are derived by applying the divergence theorem to appropriate expression on V. In 
particular see Smythe, static and dynamic electricity, page 53 eqn. 3.06(2) for a derivation of 
Green's theorem, which is, for scalars A, B and E, 



| [AV • (EVB) - BV • (EVA)]iV = J E Ay - By\s n . 8.11b 

t/ g |_ \J/l yJii _j 



Now let the function G(R,R') satisfy the equation L*G = |lR'8(R-R') in S^, where G is considered 
a function of R at fixed R'. No boundary conditions mean that G is specified to within a constant. 
Then we obtain Green's third identity: 

wW= -j G j tdSt+ jjJ^f n - G ^) dl 8.12 



80 



Magnetic fields and tokamak plasmas 



Alan Wootton 



contour I 




vacuum region 2 
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Figure 8.2. The boundaries between the plasma (S^piasma) region 1, vacuum 
(S^vacuum) region 2 and coil (S^ coi i) region 3. 

Suppose Scj) can be split up into three regions, S^piasma* S^cmim and Spoils, as shown in Figure 
8.2. Assume |l = \Xq in the plasma and vacuum region. The exterior region (the complement of 
in the right half plane) is called S^ext- Then if this external region has only linear magnetic 
material, we can apply the last equation on the region S^+S^ext- Choosing the Greens function 
Gq so that Gq(R,R') = as IRI goes to infinity, and as R goes to 0, we have 



y(R')=- \gj/s^ 



8.13 



i.e. Go(R,R') equals the flux at R' caused by a negative current at position R. Therefore we define 
¥ int (R')=-\GJ,dS, 8.14 



8.15 



and \|/ = \|/ ex t + Vint from Equation 8.12. We understand \|/ ex t as the part of the flux caused by 
currents in the exterior region, and \|/i n t is that part of the flux associated due to currents in S^. 
\|/int is homogeneous in the exterior region, and \|/ ext is homogeneous in the interior region. 



An analytic expression for Gq if |i = (lo everywhere is: 
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G (fl,fl') = — JRR' 
kit 



8.16 



where 



k 2 = - 



4RR' 



((R+R'f + iz-z'f) 

Ideal MHD 



8.17 



Here we want to note only one important equation, which is a generalization of the "virial" 
equation. We use Equations 6.1, 1.1 and 1.2. Allow the total equilibrium field to be split up into 
two parts, 



B = B 1 + B 2 



8.18 



with VxB2 = 0. We can choose the partitioning of B in a number of ways. Multiplying jxB =Vp 
by an arbitrary vector Q, we can obtain 



f( p+ jL\.Q- 



B t « VQ'Bjl 



Mo 



f 



= jT „ + JL ) (Q . n) - <»■ • Q) • <»■ ' °> k - J Q . (j x B iK V 

2jU Q J ,11,, J y 



8.19 



Mo 



with n the normal to the surface S n . We have made use of the vector identity 



Q[VxB«B] = V« 



(Q'B)B-yQ 



+ -yV.Q-B(B.V)Q 



We shall use equation 8.19 later to derive important integral relationships. 

Boundary conditions 

Last in this section we turn to boundary conditions. Suppose we have our three regions, as in 
Figure 8.2. Region 1 (S^piasma) contains all the plasma current. Region 2 (S^vacuum) contains no 
source, but contains a contour 1 on which we make measurements. Region 3 (Spoils) is outside 
region 2, extends to infinity, and contains all external currents. To find the plasma boundary we 
have to know \|/(R,z) in region 2 between the contour on which parameters are measured, and the 
plasma boundary itself. In principle we can do this knowing either 
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a) dy/dn (i.e. the tangential field) and either d\\f/d% (the normal field) or \\f on part of 1 (note 
specifying d\\f/dx is equivalent to specifying \\f to within an unimportant constant after integration. 
This is the Cauchy condition. However, there are significant problems with stability to small 
errors in the measurements. Therefore another useful set of boundary conditions is 

b) all currents in region 3, and either d\|//dn, d\|//dx or \|/ on 1. For example, suppose we know the 
currents and \\f on 1. The total fields are then the superposition of the contributions from the 
external currents in region 3 and the plasma current in region 1 : 

1//{R,Z)= W P U, Sma + Vernal 8.20 

^external is already specified by specifying the external currents, so we must only determine 
Vplasma- Since V|/ p i a sma is homogeneous outside 1, it must be completely specified by the 
condition d\\fp\ asma /dn = at infinity and either \|/pi aS ma (Dirichlet) or d\\fpi asma /dn on 1 
(Neumann). But condition 2 tells us we already have one of these specified, so that V|/ p i a sma must 
be determined everywhere outside and on 1. 

In fact we usually use an apparently over determined problem, for example knowing the external 
currents, d\|//dn and on 1. In fact this is not over determined because we only have the 

fields at discrete points, and usually the boundary conditions are only applied in a least squares 
sense. 
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9. MOMENTS OF THE TOROIDAL CURRENT DENSITY 

Let be an arbitrary function satisfying the homogeneous equilibrium equation L*0 = in S^, 
and let \|/ be the poloidal flux function which satisfies L*\|/ = -|lRj(|). Apply Green's second 
identity for the operator L* (Equation 8.11) to the pair (0,\|/):we obtain a fundamental integral 
equation: 

The moments of the current density (i.e. the integral on the LHS) are expressed in terms of \\f and 
d\\f/dn on the boundary. From L*0 = we have j 1 R" 1 (i" 1 (50/6n)dl = so there is no dependence 

upon the choice of the arbitrary constant in \\f. We introduce together with a conjugate function 
£ according to the equation 

v m = vexv^ 

VjuRJ ju 



which admits a solution L*0 = 0. In cylindrical geometry, Equation 9.2 is 



_3_' 
dR 



d 



juR) juR dz 
' 1 30 



9.3 



-, =- — 9.4 

dz\j^RJ juRdR 

The function ^ satisfies L(R- 1 (i" 1 ^) = 0, (remember the operator L = (i _1 V ((iV\)/), which reduces 
to the Laplacian if (I is uniform). Therefore the definition of ^ implies through Equations 9.3 and 
9.4 that -diR'^^l'^/dx = R^iX'^d&Zdn, where dldx is the partial derivative along 1 (clockwise on 
the outer boundary). Now by partial integration we can eliminate \\f from Equation 9.1 and write 
in terms of d\\r/dv. 

i.e. the "moments" over the current density can be measured as line integrals of the normal (B n = 
(l/R)d\|//dx) and tangential fields (B T = -(l/R)5\|//5n) along the contour 1. 
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An alternate derivation follows by letting fields q and g satisfy the equation — (V x q) = V 

ju 

so that V x (V x q) = . Then 



.r 



Jq.jrfV=Jq. Vx- dV=iV. 

v v v r 1 ' vL 



B ^ B 



1 



— xq,+— «(Vxq) dV 

\ju J ju J 



•r 



B ^ 



x qj • n + 



9.6 



This has not invoked axisymmetry. Letting q = 0V(|) = and g = R _1 £, , and using Equations 8.1 
and 8.2 gives the previous result (Equation 9.5). 

To get the notation used by Shafranov, we let q = fV(j) (which has a component only in the f 
direction, q = f/R) and consider uniform permeability (io- Then we can write an expression for 
the "multipole moment" Y m of the toroidal current density 



Y m =y\fJ,dS, =-±-j(f m B T +Rg m B n )dl 

1 iff 
where, from Equations 9.3 and 9.4, (i.e. — (V x q) = — i V x — e d 

U U K R 



= V 



various solutions of 
dR Rdz 



9.7 



—\) f m and g m are 

\juJ 



9.8 



dz RdR 

Remember that Vx(Vxq) = 0, with q = fV(|) = e^f/R, so the equation for f is 
dR 2 R dR dz 2 



9.9 



9.10 



That is, f is a solution of the homogeneous equilibrium equation and g is a solution of Laplace's 
equation. Of course, the trick is to find useful expressions for f m (equivalent to m ) and g m . 



An important point about the method of multipole moments is that the results obtained are 
sensitive only to currents flowing within the contour 1 (including vacuum vessel currents if the 
measurements are made outside this). Thus either the total equilibrium fields, or just the plasma 
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fields, can be used. The plasma fields can be calculated if external conductor currents are known. 
Using just the plasma fields alone may have advantages in terms of requiring fewer moments to 
accurately describe the data. 
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10. PLASMA POSITION 

Position by multipole moments 

Let us apply the concept of multipole moments as described in section 9 to the geometry shown 
in Figure 10.1, using the notation of Shafranov, Equations 9.7 through 9.10. Assume contains 
only the plasma and vacuum regions so that (I = (lo an d L* = A*. The term jdS^ in Equation 9.7 
is over the plasma surface, in a coordinate system centered in the plasma, while the term Jdl is 
evaluated in a laboratory frame. The two coordinate systems are related by plasma position; 
therefore Equation 9.7 must provide information on the plasma position. 




► R 




► R 



Figure 10.1. The geometry used in section 10. 

There is a contour 1, centered on R = Ri, with a local rectalinear system (t|,£), and a second 
coordinate system based within the plasma, (z,x), centered at R = R c , the plasma current center, 
defined later. Note that £ = x +Ar, T] = z + A z , R c = Ri + Ar, and R = R c + x = R\ + To be 
convenient, the moments (i.e. the choice of f m ) should, in the large aspect ratio limit 
(x 2 +iz 2 )/R p 2 «l, reduce to a form of the kind 



{x + iz) 



10.1 



with some denominator chosen to get convenient dimensions. We look for solutions of f which 
are both even (symmetric) in z, and odd (asymmetric) in z. Those suggested up to m = 2 are 
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symmetric 

/o = l 
So=0 



fl=X 



X 



1+- 

v 2R 



r 



\ 2 



1 + — 

v 2R, 



Si 



2xz 
R 



r 



x 



1 + 

V 2R c y 



( V 

2 1+— 

2z 3 



3/? 2 



asymmetric 
/o = 



/i = Z 



/ 2 = 



# 2 



2xz 



1+— 

V 2tfy 



+ ■ 



f 



V 



1+ 

v 2R cJ 



K 

2 

z_ 
R 



1+^ 



V 



Ax 2x 2 ^ 



1+ — + 



10.2 



3/? 3 



Because the Y m are sensitive only to currents flowing within the contour 1 (including vacuum 
vessel currents if the measurements are made outside this), either the total equilibrium fields, or 
just the plasma fields, can be used. The plasma fields can be calculated if external conductor 
currents are known. Using just the plasma fields alone may have advantages in terms of requiring 
fewer moments to accurately describe the data. 

We define the current center by setting Yi = 0. Using the symmetric set we obtain 



1 



2/? 



R+t 



r. 



A B 

z n 



dl = 10.3 



i.e. the current channel displacement with respect to the center of the contour 1 is 



A R = A S0 +A R1 



'so 



2R, 



10.4 



i 



10.5 



k2R, 



■ 2 r£ 
■B„+ — B. 



dl 



10.6 



Ignoring the term Aro 2 /(2Ri), then Ar is constructed from two integrals, namely 
j[^+(^ 2 /2Ri)]B T dl and j[r|+(r|^/Ri)]B n dl. The first integral is measured with a modified Rogowski 
coil whose winding density times cross sectional area varies as £, +^ 2 /(2Ri). The second integral 
is measured with a saddle coil whose width varies as T|+r|^/Ri. Alternatively the integrals can be 
constructed from discrete local measurements. 
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Consider as an example a circular contour of radius a\ based on R = Ri Then \ = a\ cos(u)) and T| 
= ai sin(Q)). These coils measure the fields in the coordinate system (p,co,<f)) based on the vessel 
center, and B n = B p , B T = B ro . Assuming a plasma with no vertical displacement, symmetric 
about z = 0, we can write 

B m = B T = + cos(n(o) ] 10.7 

2m, L „ J 

B p = B n =^Wju nS m(na>)\ 10.8 

Substituting these expressions into Equations 10.4 to 10.6 gives 

A R ~^(A l + J u l ) + Ml + ^ + J u 2 )--^-(A l+J u l ) 2 10.9 



4ft, 



8ft 



i.e. for small displacements, (k\+\i\) « 1, and a circular plasma (X n , (l n , = for n >1) we have, 

A R ^^ L (A l +jU l )+— 10.10 
r 2 v i f-u 4 ^ 

i.e. all that is needed to measure the displacement of a nearly circular plasma within a circular 
contour (k n , |i n , n > 2 = 0) is a modified Rogowski coil whose winding density varies as cos((fl), 

and a saddle coil whose width varies as sin((0). This simple coil set gives the correct answer 
when the constant ai 2 /(4Ri) is allowed for. We already knew this. To allow for significant non- 
circularity the more general expression (Equations 10.4, 10.5 and 10.6) should be used. We can 
also derive equations to determine the vertical displacement of an arbitrarily shaped plasma, 
using the asymmetric components in Equation 10.2. 

Application to the large aspect ratio circular tokamak 

Let us apply these ideas to the circular equilibrium described in section 6, surrounded by a 
circular contour on which we have a sinusoidal area Rogowski coil and a cosinusoidal width 
saddle coil. The equations given in section 6 were in a coordinate system (r,9,<f)) based on the 
plasma geometric center; they were transformed into the vacuum vessel coordinate system in 
section 7, Equations 7.1, 7.2 and 7.3, allowing for a geometric shift A g . 

The output from the integrated saddle coil with n w layers of width w((fl) = wosin(G)) and integrator 
time constant Tint is 
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7 int ^ 



int 



Using the expression for B p from Equation 7.3 we get 



a, 



r 



( \ 



a, 2R, 



In 



+ (a4 



1 1 



V a, 



10.11 



10.12 



The output from the integrated 'modified Rogowski coil each turn of area A, with nocos(O)) turns 
per unit length, is. 



i = ^|\ ( ^ )cos( ^ = iv^ /l 



int 



10.13 



Using the expression for B ro from Equation 7.2 we obtain 



3 



\ 


a, 






r 


( 



r r ^ 



In 



A+ln 



1 



+ IA + - 



10.14 



Before we can substitute these expressions (Equations 10.13 and 10.14) into Equation 10.10, we 
must recognize that our equilibrium fields were evaluated in a left handed coordinate system, 
while this section we have worked in a right handed system. Sorting this out we find A,i => 
and (ii => Hi, so that 



a\ ( 1 



2R, 



10.15 



This is the difference between the geometric center A g and the current center Ar of a circular 
plasma, under the present approximations. We also note that, after sorting out the coordinates, 
subtracting the outputs from our coils gives 



2 // --2. 





f \ 




A + ln 




V 




) 



10.16 



that is, we can measure 
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#+-^+1=1+^-/0-111 

2 a, 



f \ 
a 



10.17 



Combining Equations 10.15 and 10.16, we see that we can measure the current center, the 
geometric center, and A, using our modified Rogowski and saddle coils. 
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11. PLASMA SHAPE 

Using higher order moments we can obtain information on the plasma shape. Y2 determines 
ellipticity and Y3 determines triangularity. Using equations 9.7 and 10.2, we obtain: 



Y 2 



V 



R 



1 J 



^ A ^ 

v R\j 



11.1 



where 



1+ — 

R x j 

r 



R^JA 



TL V 2R X 3R& 
= jB T dl =ju l p 



B T d\ 



11.2 

11.3 
11.4 



That is, with the Rogowski coil measuring I p (i.e. sq 1 ) and either modified Rogowski and saddle 
coils, or single point measurements of B n and B x suitably combined, we can construct Y2. If we 
want to use modified Rogowski and saddle coils, then to obtain I p , Ar and Y2 takes a total of 5 
coils. For a circular contour, and ignoring toroidal effects, Equation 11.1 is written as 



Y 2 =-A 2 R +A 2 z +^(Z 2 +jU 2 ) 



11.5 



That is, neglecting toroidal effects we need only and (12, in addition to Ar and I p . 

To interpret the moments it is necessary to assume a plasma current distribution; because the 
moment is an integral of the current density over the surface there is no unique solution for the 

boundary shape. As an example, consider a uniform current density and a surface described by an 

(A 2 (z\ 2 

ellipse with minor and major half width and half height a and b, so that 1 — + 1 — = 1 . Then 
b 

for k = 1, k « 1, and ignoring toroidal effects, we have 

a 



ka 



11.6 



In a similar manner, if the surface is described by a function which includes both elongation and 
triangularity, namely 
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+ 



a 



2 '\ 



& ) 



11.7 



then ignoring toroidicity (i.e. x/Rj and z/Ri «1, we have for a flat current distribution 



11.8 
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12. MOMENTS OF PLASMA PRESSURE 

The Virial Equation 

Here we consider certain integral relations which allow us to determine information on the energy 
density within S^. The techniques described to determine plasma position, and more generally 

plasma shape, needed only Maxwell's equations. Now we add the constraint of equilibrium. 
Generally we apply these results to = S^piasma + S^vacuunv even simpler relations are obtained 
by identifying with S^pi aS ma- 

From a previously derived integral equation (Equation 8.19), we can obtain some very useful 
relationships. First, take B2 = (i.e. all the fields are described by B = B1), and let the arbitrary 
vector Q be the polar vector (r = ReR + zez-)- Then Equation 8.19 becomes the "virial theorem": 

ff B 2 ) M 5 2 V x (i-B)(B^n)l 

\3 P+ —dV= \p + —Uvn)-± * ys n 12.1 

Note that if B were the self plasma field, and there was no other field, then the surface integral on 
the RHS approaches zero as the surface (where p = 0) approaches infinity, because B(self) R~ 3 . 
This contradicts the positive definiteness of the LHS volume integral, showing that equilibrium 
by self fields alone is impossible. 

Now restricting to toroidal symmetry, d/d§ = 0. Therefore the poloidal (subscript p) and toroidal 
(subscript (])) fields satisfy 

B = B p + B, 

„ 12.2 

Next assume that the surface S n is outside the plasma, where VxB = 0, and that the vectors Q and 
B2 are purely poloidal (i.e. there is no toroidal component of B2); 

Q=Q R e R + Q z e z ^ 

B 2 = B 2R e R + B 2z Q z 

Noting that the toroidal field outside the plasma is given by B^e = const/R, then we have 



J < (Q • n)dS H = J<V • Q - dV 12.4 
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and 



B 1 *B 1 .VQ=B pl .B pl .VQ+ ^ 



12.5 



Using these, we can express Equation 8.19 so that the toroidal field enters only as B^-B^e 2 on 
the LHS: 



R 2 A 



Mo 

(B. •Q¥B.«n)l rr , xn 

A) 



12.6 



Now let 



T = p + 



2jU Q 



_£l is. i 



v 2ju 







12.7 



12.8 



(Remember Bi is the poloidal field). Then the integrand on the LHS of Equation 12.6 can be 
written as (from now on drop the subscript p for poloidal) 



L '" J IdR dzJ R^ R 2 Mo Vdz dR) Mo *<dz BR) 



12.9 



Now we have to chose something for Q. Let 
u = R- R ch +iz = pe im 



12.10 



with R c h some characteristic radius, a fixed point within the plasma cross section. Then 
p 2 =(R-Rj+z 2 12.11 



Now with this notation let 

Q R = F(u) 
Q z =-iF{u) 



12.12 
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so that the last two terms in Equation 12.9 are zero. Next take 



F{u)= u 



12.13 



Finally we can rewrite Equation 12.6 in the form 
\^2(m + l)Pu m +^u m+1 ^iV 

= f[p + £] (n s - Mj- ( ^-^ )(B ' ,n) K - jp • (B 2 x Q)> 

4"(K Mo v 



V 



12.14 

Surprisingly, this is useful. We now consider the results when we take different values of m. 

m = -1, B 2 = B ext 

Take B2 as the externally supplied maintaining field and Bi the plasma self field. For large 
volumes the surface integral approached zero, and we have 



12.15 



or equivalently 



\j,B z>ext RdS,=-j 



p + 



' B 2 pl+ Bl-B lM 



v 



2ju J 



12.16 



Suppose (low toroidicity) that p/R c h «1. Then the R dependence of T appears in the form R - 
R ch , so that dT/dR = -dT/dR ch . Now because T/R = V Te R - dT/dR (if p/R ch «1) we have, with 
T = at infinity: 



{L dv= (±L dv «([ Tdv ) 
[r I dR dR ch V > 



12.17 



Now introduce some definitions: 



B, = 



2nR 



12.18 



L cA p V 



96 



Magnetic fields and tokamak plasmas 



Alan Wootton 



\TRdS, 

R T =— r 12.19 

\TdS $ 



L=—A^dV 12.20 



pt 



fr=-^j\pdS^ 12.21 



Mi 



dS< 12.22 



Note that (ioLpt/(47t) is a total inductance, and the usual "inductance per unit length times 47l/|io" 
is given as 

,2 



L = A^dV 12.23 

p u n Rjt J 2u n 



Now we can write Equation 12.16 as 

This is a generalization of the external vertical field needed to maintain a plasma, allowing for 
non circularity. 

m = 0, B 2 = 

Now Equation 12.14 gives an integral expression along the minor radius: 
(Rr-R ch )i 



\PdV+ KT ch, [TdV 



V T V 



= 0.5 \ I V \ " lJ n • -2*2* t • e p ^S„ - 0.5 J j, [* z2 (/? -R ch )- B R2 z}lV 



12.25 
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If we assume that IR-R c hl « R c h, then the term with (Rx-R c h)/(2RT) can be omitted, and the LHS 
is independent of the choice of S n . Therefore the RHS must also be independent of S n . Then we 
can write (with IR-R c hl « R c h) 

j Pdv =^hll Si 1226 

V ^ 

where si must be independent of the choice of B2. In particular, taking B2 = 0, Equation 12.25 
can be written, using the definitions above, as 

/3 I =ju I + s l 12.27 

where 

= \ IK - ^ > n ' e , - 2B r B JK ' 12.28 

The difference between si and 1 is due to non circularity. 

In fact, the constraint R-R c hl « Rch, but not (RT-Rch)/(2Rx), can be relaxed by redefining |3r v 
and (j-iv as volume integrals: see Equations 12.31 and 12.32 below. 

m = -1, B 2 = 

Now we obtain another important integral relationship, assuming IR-R c hl « R c h, namely 

R + h. = h + h 12.29 

Hl 2 2 2 



where 



=lh \ - 5 » > • e , - 2B ^ • % } S n 12-30 



and C, = R - R c h. 

In fact, the approximation IR-R c hl « Rch is not necessary if (3 and (J, are redefined as volume 
integrals, namely 

27tR c J IV =^- 2 \pdV 12.31 
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8x f - £jL) 
2ju 



2xR chJ u IV = -- "^J " ^ dV 12.32 
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13. Pi + li/2 
Solution 

In the previous section, Equation 12.29, we showed that we could express the parameter Pi + Lp/2 
through the integral equation 

j 3+ Ll = ?i + s 131 

Hl 2 2 2 

where si and S2 are two integrals of the fields over dS n , which we will be able to measure as 
contour integrals, namely 



' ' ~^br J [( B l ~B 2 n )pn'e p - 2B T Bjn • e p ]lS n 

R I,' „ r- 



w;* l32 

^ K - f [fe 2 - 5„ 2 Vn • e„ - 25^ jn • e „ "fcfl 



S 2 ..2 



**ch*p 1 



13.3 

Equation 13.1 was good for toroidal geometry if volume definitions (3iv and (iiv were used. Note 
that dS n = 27lRdl, so that si and S2 can be written as contour integrals around 1. Unfortunately 

they involve the squares of fields, and so we cannot design simple modified Rogowski and saddle 
coils to make the measurements. Instead we must measure B n and B x at discrete points along the 
contour, and then construct the required integrals. All we have to do is construct the integrals si 
and S2 from measured B n and B x . 

Suppose we are measuring fields on a circular contour of radius ai, centered at R = R\ = R c h; that 

is we identify our characteristic radius with the center of the contour 1. Then in the vessel 
coordinates we have x-e p = 0, n-e p = 1, n-e^ = cos((o), x-e^ = -sin(co), and dS n = 
27lRi[l+(ai/Ri)cos(G))]aidG). If the fields can be expanded as in Equations 10.7 and 10.8, and 
keeping only the terms A,i and \x\ (i.e. a circular plasma), we obtain 
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s 2 =^(A l +ju 1 ) 

a. 



13.4 



i.e., substituting 13.4 into 13.1: 

/? /+ -^ = l+% 1+A ) 

2 a, 



13.5 



Let us apply this equation to the displaced analytic equilibrium discussed in sections 6 and 9. 
The Rogowski coil with cosinusoidal varying winding density, and the saddle coil with sinusoidal 
varying width, tell us Xi and First there is a mess with right and left handed coordinate 

systems to unravel. Doing this then (I => -|i in Equation 13.5. Now Equation. 10.16 already tells 
us that 



r 



In 



Hl 2 



1 



13.6 



Therefore comparing Equations 13.5 and 13.6 we must have (allowing (I 
that 



-|l in Equation 13.5) 



l± + In 

2 



\<* P J 



2 
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This is exactly what we expect, because the total inductance to a radius ai is given by: 



( Mo 



= UaR 



r 



In 



a, 



2 
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Separation of ft and I, 



If we have Pi + 1/2 from the poloidal field measurements just described, and Pi from diamagnetic 

measurements (see later) then obviously we can separate Pi and \\ If no diamagnetic 

measurement is available, two possibilities exist. For non circular plasmas, there is a third 

integral relationship, which I have not derived, which gives in terms of a measurable line integral 
the parameter l v (2p+B z 2 /(io)dV. When V is the plasma volume, this is related to 2Pi + L p . If the 

volume averages <B Z 2 > and <B p 2> are different, as is the case for non circular discharges, then 
this measurement allows the separation. 
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For near circular plasmas, we must estimate Lp separately. For example, for the simple circular 
low beta equilibrium of section 6, we can take a model current distribution jtyo(i) = jo(l-(r/a) 2 ) a . 
Then lj = Lp-ln(ai/a p ), with h given as a function of a = (q a /qo -!)• By assuming qo = 1 we can 
then estimate lj, and make the separation. 

Comments on the definition of poloidal beta 

We must be careful with the definition of "poloidal beta". So far we have used Equation 12.21, 
namely 



We could replace (3i by (3 p (the poloidal beta), which characterizes the ratio of plasma pressure to 
the pressure of the magnetic field for an arbitrary shaped cross section. It should be introduced so 
that the pressure balance, Equation 12.27, is replaced by 




13.9 




13.10 
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14. DIAMAGNETISM 

Comments 

Next we turn to the toroidal flux and its measurement. Before plasma initiation this is simply 
given by ^B^dl = Mo\ J z dS , i.e. IkRB^ =jU q I z , where we have applied Amperes law to a contour 

in the toroidal direction, encircling the inner vertical legs of the toroidal field coils (i.e. I z = 
i z n t n c , with i z the current from the generator, n c the number of assumed series coils, and n t the 
assumed series number of turns in each coil). 

Microscopic picture for a square profile piasma in a cylinder 

During formation inside a magnetic field the plasma particles acquire a magnetic moment m: 



m = Area orhiI I orhiI = 4 — ] | — ] 14. 1 



Since ft) = qB/m e , we have 

m v 2 

m = ^— 14.2 
IB 

adding up to a total magnetic moment 

M=nmS lp 14.3 

per unit length of column with cross section and a number density of n Supposing cylindrical 

geometry the elementary currents cancel within the homogeneous column, leaving only an 
azimuthal surface "magnetization' current density j s : 

nk b T p, 

J= nm = — h — = J ^ L 14.4 
B B 

where pj^ = nkb(T e + Tj)j_, kb is Boltzmann's constant. The toroidal field will be modified. 
The associated flux from this surface current can be calculated: 
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A<D = ^a 2 AB s = m 2 -^-I = xa 2 n -&-2xRj, 
" * p 2nR p 2kR 

p 2nR B a p ^° B a 



Using the definition of Pi, discussed much more a little later, we then have 
A4> = ^° pHl 



Macroscopic picture 

Let us consider a toroidal device with no toroidal current plasma current, i.e. a stellarator, in 
which the necessary rotational transform is produced only by external conductors. Starting with 
the radial pressure balance, with p^= at the plasma edge, and approximating the torus by a long 
cylinder, then 

^ = JA 14-5 

integrating over the minor radius (r = to a p ) gives 

a ,> 

P L = -B,\j e (r')dr< 14.6 



and 



^ a ,> ft a P °p 
(p L ) = — j \2np 1 rdr = \27trdr\ j (r')dr' 

7777 •> 7777 J J 



14.7 



B "" a "S(h 
= M m- 2 j s (r)dr = -5 \-^j s (r)dr 

. ap 

Then j se = J S(r)/(7ta p 2 )je(r)dr is the effective surface current density at the plasma edge as a 
consequence of the finite plasma pressure. 
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Paramagnetic and diamagnetic flux 



Outside the plasma the toroidal field has the form B^e = B^Ro/R), with B^q the value at a fixed 
radius Rq. This toroidal field, together with the poloidal field, takes part in balancing the plasma 
pressure. 

We need an equation relating fields to pressure. Substituting V x B = ju j into Vp = j x B yields 



(using Bx(VxB) = 



VB 2 



(B • V)B ) 



B 2 



p + 



(B-V)^- 



For a straight axially symmetric system (d/dz = 0) we obtain 

f ->2\ 



dr 



B 2 + B e 
- 

2ju J 



v 



JL 

Wo 



Multiplying each side by r 2 , letting u = r 2 , du = 2rdr, dv = d/dr(.), v = (..), we obtain by 
integrating by parts (Judv = uv - jvdu) 



Bl +Bl) 
2ju J 



Bl + B: 



D T D 

P + 



B 2 

2rdr = - \—rdr 



i.e., 



B 2 + Bl 
9 - 



V 



2ju 



= — f 



B 2 ) 
p + — YlTCrdr 



That is, ignoring curvature and equating B z with B^, the pressure balance constraint is 

2ju ( P )=K+K-( B D 148 

where B^ is the toroidal field inside the plasma, B^e is the toroidal field outside the plasma, <..> 

means an average over the plasma radius, and we have assumed p = at the boundary (i.e. at r = 
a). That is, for a given plasma current I p and pressure <p>, the difference (B^e 2 - <B§ 2 >) adjusts 

itself to ensure pressure balance. This happens because of a poloidally flowing current, either 
diamagnetic or paramagnetic, in the plasma, as we derived in equations 14.4 and 14.17. In a 
tokamak we have B^e 2 » B^e 2 - <B ( ) ) 2 >, so that if the cross section is circular with radius a p , 
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where 



<-(^ 2 )-25^-^) = ^p 14.9 



S^=m 2 /\B^-B^ 14.10 



is the diamagnetic flux of the longitudinal (toroidal) field. We will discus its measurement later; 
it is the difference in toroidal flux in the plasma column when the plasma is initiated. Defining 



i.e. B, = f or a circular cross section 14. 1 1 

B 

da 

with B 8a = ^° p the poloidal field at the plasma edge, we can write 
SxB.S® 

A = i+— £5- 14-12 

From this equation we write the net flux difference 8$ = (|loIp) 2 /(87i:B ( |) e ).((3i - 1) as the sum of 
the paramagnetic flux 84> p : 

2 2 

P 8^ 

due to the poloidal component of the force free plasma current, and the diamagnetic flux 83>d: 
d^ d = -d^J3, 14.14 

due to the poloidal currents providing pressure balance for the finite pressure. 

Toroidal, non circular geometry 

In a torus curvature must be accounted for: Corrections with coefficients (a/R) appear in the RHS 
of the equation for Pi. For (3i « R/a the corrections are small, and it should be noted only that 
for B^e the value on a line R = R g , the plasma center, should be used. We will derive that below. 
Also, we must now be rather careful with definitions. 
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We derived, in section 13, a general expression relating Pi to poloidal field measurements and the 
diamagnetic measurements. This derivation allowed for both toroidicity and non circularity, and 
gave in terms of volume integrals (equations. 12.31, 12.32, 12.27, 12.28) 

J3 iv =ju i +s 1 14.15 
M„ = A —\&^dV 14.16 

& v= vihi pdv 14 - 17 

The measurable integral si = 1 for a circular discharge. Remember that R c h was a characteristic 
major radius of the plasma. Because we can write IB^-B^I « B^o, the diamagnetic parameter 
can be expressed through the experimentally determined flux (the measurement is discussed later) 

S*=\(B+-B,)dS, 14.18 

By substituting (B^-B^ 2 ) ~ 2B^(B^ e -B ( ) ) ), and writing 

B Ul9 

* R 

then we have 

8xB M 5<f> 

M 1V = f- 2 14-20 

That is, in toroidal geometry and for an arbitrary shaped cross section, measuring 8$ and the 
integral si gives us a good measure of Piv, the volume defined current beta. Note that the field 
B^o is the external field at the characteristic plasma major radius, which we can take as R c h ~ R g .. 

The meaning of pi 

We can understand Pi by integrating by parts Equation 14.1 1 to give (for a circular cross section), 
and assuming p at r = a p = 0, we have 
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ItPp 



AVo o o AVp o 



r, 



\pr 



2\"p 



[r^dr 
J dr 



p o 



14.21 



8^ Lidp 
Moll ' 



Here we have integrated by parts, with u = p, du = (dp/dr)dr, v = r 2 , dv = 2rdr. Next we use the 
approximate pressure balance equation (dp/dr = jeB^ - j^Be, we will work in a RH coordinate 
system). We substitute for je and in terms of and Bq (from ^ j=VxB), i.e. 
dB A 1 d(rB g ) 

dr 



MoJe = 



or r 



, we obtain 



dp d 
dr dr 



B 



+ ■ 



B a d 



\2jU Q J jU Q r dr 



K)=o 



14.22 



Substituting for dp/dr from equation. 14.22 into equation. 14.21 gives 




14.23 



If (dB^ 2 /dr) > then pi > 1, and if (dB^/dr) < then pi < 1. Figure 14.1 shows this 
schematically. A reduced B^ (diamagnetic) inside the plasma is associated with Pi > 1. The 
normal ohmic heating case, with Pi ~ 0.3 has an increased B^ inside the plasma (paramagnetic). 
The figures are not well drawn, because we must have that B^R is constant on a flux surface. 
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B 



P, >1 



B 




R R 

Figure 14.1. The toroidal field with different values of (3t. 



Measurements 




coil output 

Figure 14.2. A simple loop to measure toroidal flux. 

To measure 84> it is usual to use a wire wrapped around the vacuum vessel, as in Figure 14.2. 
The voltage output from the coil is integrated to give a signal proportional to the enclosed flux. 
For most cases 8$ is small, typically ~ lmWb, which must be compared to a typical vacuum flux 
<£ v = lWb enclosed by the same loop. Therefore we have to use techniques which allow 
measurements to better than 1 part in 10 4 to get 10% accuracy in values of Pi. This is done using 
compensating coils. 
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loop 2 



Fig 14.3. Diamagnetic loop Fig 14.4. Two concentric diamagnetic loops 

and compensating coil 

Figure 14.3 shows a typical set up for the 'diamagnetic loop', a single turn coil around the vessel, 
and a compensating coil. The idea is to make the vacuum signal from the two coils as near 
identical as possible, so that a simple summing (or subtraction) circuit can be employed. A 
compensating and balancing circuit is then employed, as for example is shown in Figure 14.5. 
The vacuum common mode flux is canceled using the summing integrator, so that the balance is 
determined only by the effective areas and the resistors. The product of resistance and 
capacitance (RC) is adjusted using 'toroidal field only' shots (no plasma). By opening the switch 
after the toroidal field (TF) has 'flat topped' much of the common mode is also rejected. Of 
course, to balance the system during the toroidal field only shots the switch is opened to consider 
toroidal field ramp up. The circuit shown also allows for phase difference between the two 
signals to be compensated. 
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plasma 




diamagnetic loop 



compensating coil 
I 



balance resistor 



integrator gate 




Figure 14.5. A diamagnetic coil compensating system 

An alternative technique is to wrap two simple toroidal flux loops around the vessel, but at 
different minor radii. Such a configuration is shown in Figure 14.4. Let two concentric loops 
have radii bi and b2, and let R\ be the major radius of the loops. Then, remembering that B^e = 
B^oRi/R, we have after time integration 



14.24 



<!>(&,.) = 2;^ 



R 1 -(Rf-bff 



14.25 



Si>=<t>(b l )-k(<i>(b 2 )-<t>(b l )), 



k = 



b 2 JW^+Jr^ 



14.26 
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From these two signals we can calculate 8$. The subtractions are performed electronically (i.e. 
analog); the constant k is determined experimentally so that 8$ is zero without plasma (toroidal 
field only). 

In reality there are problems. One in particular is caused by the discreteness of the toroidal field 
system, and the current redistribution in the toroidal field coils during a shot. When the toroidal 
field current is initiated, the current flows at the inner edge of the conductors to minimize the 
linked flux. As the pulse proceeds, the current redistributes and approximates a uniform 
distribution (not exactly because of repulsion of current channels). The time for this 
redistribution to occur is approximately the radial penetration time of the poloidal current into the 
conductor of radial extent w: x ~ 7t|lo<3w 2 /16, typically 200 ms. If the toroidal field system was a 

perfect toroidal solenoid this redistribution would leave the fields unaffected. However, because 
of the discrete number of toroidal fields, there is now a time varying toroidal field ripple. The 
size of the changing field ripple depends on where a pickup coil is placed: therefore two coils 
linking the same steady state flux can link different transient fluxes. It is best to place the coils 
between the toroidal field coils, where the redistribution effect is smallest. 

Another problem is due to poloidal eddy currents in any conducting vacuum vessel. This 
produces a non zero change in the toroidally averaged B^, not just a local ripple. Therefore it 

couples strongly to the pickup loops. Compensation for both the effects discussed has been 
performed successfully using software, by simulating redistribution and eddy currents as simple 
circuits, coupled to both the primary coil current and a secondary pickup coil. 

Further problems occur if the loops are not exactly positioned, so that they couple to the poloidal 
fields produced by the primary, vertical field and shaping windings. That any such effects exist 
can be checked for by firing discharges with positive (+) and negative (-) B^. Let A<£ p be the 
signal caused by poloidal field coupling, A<£(+) the signal obtained with positive B^, and A<£(-) 
the signal obtained with negative B^. While the toroidal field coupling effects will change sign 
with reversing B^, any poloidal field effects will not. Therefore 

A<D(+)-A<D(-) = 2((53> + A<D p ) 14.27 

A4>(+) + A<D(-) = 2A<D„ 14.28 

Any finite A4> p can be correct for using a circuit model again, with the coupling between any 
winding (including vacuum vessel) and the pickup coil written in terms of mutual inductances. 
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15. FULL EQUILIBRIUM RECONSTRUCTION 

The problem under discussion is how to reconstruct as much as is possible about the equilibrium 
from external measurements. In particular, if we knew everything about the fields outside the 
plasma, what could we uniquely determine? Could we separate (3i and lj? (In principle yes for 

toroidal systems). Could we go further and actually uniquely determine the plasma current 
distribution of a toroidal plasma (I don't know)? 

If we want to allow measurements of dp(\|/)/d\|/ and FdF(\|/)/d\|/ we must have redundancy beyond 
that required to solve the equilibrium equation with a known current density profile. We see 
immediately a problem in straight geometry, because with straight circular cross sections the 
magnetic measurements must be consistent with a solution that has concentric circular surfaces. 
In this case the measurements are consistent with any profile function which gives the correct 
total plasma current, so there is an infinite degeneracy. 

In practice the mathematical subtleties of what can and cannot in principle be determined are not 
discussed. Instead the usual technique for equilibrium reconstruction is to choose a 
parameterization for j^, with a restricted number of free parameters. These free parameters are 

chosen to minimize the chi squared error or cost function between some measured and computed 
parameter, for example the poloidal field component on some contour. Obviously we need some 
boundary conditions, as discussed in section 8. Information available might consist of the flux \\f 
on a contour, the fields B n and B x on a contour, and currents in conductors, the total current, and 

if lucky the diamagnetism. Typical parameterizations are 



g{§¥) 15-1 

= (a x 8y/ + a 2 8y/ 2 +a 3 Syf)R + (b l Sy/)R~ 1 15.2 

with 8\|/ = (y-yboundary)(Vmag axis"¥boundary)- m each case the term proportional to R represent 
the part of proportional to dp/d\|/, and the part of proportional to R _1 is proportional to 
FdF/d\|/. R G is some characteristic radius. The assumed function g might be of the form g(8\|/,y) 

= exp(-y 2 (l-8\|/) 2 , d\|/Y or 8\|/ + y8\|/ 2 . If in equation 15.1 the function g(8\|/) is 1, then the 
description of is called quasi uniform: for a circular outer boundary there is an exact analytic 

solution of the equilibrium equation. 
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The general results show that, given only poloidal field magnetic measurements (i.e. no 
diamagnetic signal), then (3j and lj can only be separated for significantly non circular equilibria. 
From our discussions of Pi+h/2, we know we can measure (3r+Lj/2 exactly anyway. If we have 
reconstructed the plasma surface then we can evaluate the integrals si and S2 on that plasma 
surface, and then Lj = lj. For circular equilibria it is found that only the sum (3i + \[/2 can be 
found. Of course, adding a diamagnetic measurement allows the separation for any shape, as we 
have seen. Also adding a pressure profile form, or the pressure on axis, is enough to allow 
separation even for circular discharges. It is generally found that a distortion of the shape to b/a ~ 
1.3 is needed to make the separation. 

A significant problem is that, having made specific assumptions about the form of the current 
density profile, we do not know how general our results are. 



114 



Magnetic fields and tokamak plasmas 



Alan Wootton 



16. FAST SURFACE RECONSTRUCTION 

Here we want to investigate 'fast' methods of determining the plasma boundary, and fields on that 
boundary. If we can do this, then we can exactly derive parameters such as (3j + 1/2, and q, on the 

boundary, as well as parameterizing the boundary shape itself. Here we specifically do not want 
to obtain any information on fields and fluxes inside the plasma boundary. 

In the vacuum region S^vacuum bounded by the contour lpiasma and 1, the flux function satisfies 
the homogeneous equation L*\\f = 0. A solution for this equation which agrees with boundary 
conditions on 1, and that is valid in a region containing S^vacuum* must then be valid on the 

plasma boundary itself. 

In general we choose to approximate by a series solution 

N 
7=1 

With \|/° any known terms, and the basis functions all satisfy the homogeneous equilibrium 
equation on some region S^o which includes the vacuum region S^vacuum- Following the 
discussion in section 8 we can write y = V[/ ex t + Vint from Equations 8.12 through 8.15. Usually 
we take an interior solution \|/i nt associated with current distributions inside S§, and an exterior 
part \|/ ex t associated with currents in the exterior region S^ext. For example V|/ ex t rnay be known 
from measured currents in external conductors, or calculated directly from the measurements 
using Green's functions. Various representations for the plasma current have been used; 

an expansion in toroidal eigenfunctions, 

discrete current filaments, 

single layer potential on a control surface, 

To date no purely analytic answers are available; generally the coefficients in an expansion, of 
position or currents in filaments, are altered numerically to provide a minimum "chi squared" 
between some measured and computed fields or fluxes. For example, suppose the actual 
measurements are represented by qi, with i running from 1 to M. There is a response matrix Q 

such that 
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where q es t,i is the value expected for qi when Equation 15.1 holds. q es t ,i is associated with V|/°, 
and the known Qy is associated with %j. Then the usual least squares approach will determine the 
coefficients cj to minimize the function 

v 16.3 

.■=i °7 

With Oj the standard error of the i th measurement. This procedure may not be stable, in which 
case some numerical damping is added 

One technique which avoids the iterations necessary to match the measured and predicted field 
components is as follows. From the measured fields construct the multipole moments Y m using 
the techniques outlined in sections 9 and 10. Then m toroidal filaments with current (I p /m) can 
be positioned to give the same moments Y m as those measured. Because we have analytic forms 
for the Y m produced by discrete current filaments (the integral Y m = Jj^fmdS^ only takes a finite 

value at the filament location) we can derive analytic expressions for the filament positions in 
terms of the measured Y m 's, thus avoiding the need for the iterative procedure. Just as discussed 

above, we then use toroidal filaments with known currents for external windings, the m filaments 
for the plasma, and plot the flux contours immediately. Solutions up to the plasma boundary are 
as exact as our set of moments allows. I am not sure how unique the solution is, or to what extent 
I should consider taking more than m filaments. We have still to ask if our solution is unique: 
that is, do the m moments uniquely specify the fields on the contour 1? 

An example of such a procedure is shown in Figure 16.1a. Figure 16.1b shows a full 
equilibrium reconstruction with j^r) iterated until a good fit between measured and computed 

moments was obtained. Clearly the 3 filament approximation, with the filaments chosen to give 
the measured moments Y1Y2 and Y3, gives a good description of the outer surface. 

In principle we should be able to extend the "moments with filaments" method of finding the 
plasma shape to the use of an analytic representation for the current density. Indeed, we did this 
in section 11 to find a relationship between the second moment Y2 and ellipticity. However, 
there we made an arbitrary choice that the current density be flat. In fact we should specify that 
satisfy the Grad Shafranov equation: this problem is considerably more complicated. However, if 
solved, we should be able to obtain analytic relationships between the measured moments and the 
plasma shape. We would still have to parameterize the form for j^: this would be restricted by, 

for example, a knowledge of the ratio of q on axis to q at the edge. 
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with the experimentally measured moments. The Flux surfaces predicted by a free-boundary equilibri, 

Mm line represents the plasma surface, the circle the calculation. Flux-independent current, with poloidal beta - 

"•nwm vessel around which the line internals were measured. gives moments similar to those measured 



gure 16.1a. External surface reconstruction Figure 16.1b. External and internal 

surface reconstruction 
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17. FLUCTUATING FIELDS 

(MIRNOV OSCILLATIONS and TURBULENCE) 

Mirnov Oscillations 



In tokamaks it is expected that magnetic islands play a role in determining transport. Their 
structure is approximately of the form exp(i(m0+n(|))), and they are located at surfaces where q = 
m/n, a ratio of integers. This is illustrated in Figure 17.0. 



resonance 1 
mi 




n l 



- resonance 2 



Figure 17.0. The envelopes of two adjacent magnetic islands. 



In terms of the Fourier coefficients of the radial component of perturbing field b r>mn at the 
resonant surface (where q = m/m), the island full width w is given by: 



17.1 



Note that if the current distribution is uniform, (so 5q/6r=0) and a resonance exists, this predicts 
that an infinitesimally small perturbing field will destroy the circular flux surfaces completely. 
Mirnov first studied their presence using be loops, actually measuring db$ldt outside the plasma 
(B denotes total fields, while b denotes just the oscillating part). Because the coils are outside the 
plasma they do not measure the field strength at the integer q surface where the perturbation is 
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resonant, so that we cannot immediately calculate the magnetic island width using Equation 17.1. 
Other problems to cope with include how to determine the poloidal (m) mode number from a set 
of coils measuring be if they are not on a circular contour, so that a simple Fourier analysis is not 

possible, and toroidal effects. 



If the field perturbation is b, then outside the circular plasma at some radius r > a p we have 

17.2 



„ , ob r im , in, 
V. h = -L + — b+ - b=0 

or r R 



(Vxb>e r = - b,--b g 
r R 



Form the second equation we have 
mR 



17.3 



17.4 



nr 



which for low m, n is »1, i.e. be » b^. Therefore with Equation 17.2 we determine that the 
most important contributions to measure are b r and be, not b^. If we are measuring just inside a 
conducting wall (e.g. the vessel) then b r = 0, so that only be should be monitored (but note the 
conducting wall also affects the value of be: it can also affect the instability which produces be 
itself). 



coil at r 




^^•Jilament at r = a ,0 = 0o 




Figure 17.1a 



Figure 17.1b 



A common representation of the fluctuations is as current filaments aligned along the field lines. 
Initially let us consider the fields produced in a straight cylinder by a current filament (current I) 
aligned along the cylinder, located at poloidal angle = 0q, at r = a (see Figure 17.1). Then. 



b e = 



jU Q I (r-acos(0-0 Q )) 
In (a 2 +r 2 -2arcos(6>- #„)) 



17.5 
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r 2 K (a 2 + r 1 - 2ar cos(#- O )) 

These can be rewritten by expanding the denominator into [r - a exp(i(0-0o))].[r - a.exp(- 
i(0-0o))]. For example: 

h &l\ 1 i 1 I 177 

These new expressions can be expanded in powers of (r/a), so that 

r < a r > a 

inside filament outside filament 

27tr n=l ^a J 2^r-nVr/ 



n=Q 



17.8 



i.e. we have the fields as a Fourier series. A certain sum of the coefficients allows us to look only 
at the current inside the cylinder, while ignoring external currents: 

27t 27t / CL\ m 

r^b e {0)cos{mO)dO-r^bXO)s\n{mO)dO= Mo\~J cos(m# ) if a < r (outside) 

oo r 

lit lit 



rjb (O)cos(mO)dO-rjb r (O)sin(mO)dO= if a > r (inside) 



o o 

In In 



r^b e {0)sm{mO)dO+r^b r {0)cos{mO)dO = Mo\~ ) sin(m# ) if a < r (outside) 




lie lie 



rjb d (O)sin(mO)dO+rjb r (O)cos(mO)d0=O if a > r (inside) 



17.9 



Consider just the m = 1 component . These equations tell us that, in straight geometry the sum of 
a modified Rogowski coil with winding density proportional to cos(0), and saddle coil with width 
proportional to sin(0), tells us the horizontal position. A modified Rogowski coil with winding 
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density proportional to sin(0), and saddle coil with width proportional to cos(0), tells us the 
vertical position. Actually we already new this; section 9 tells us how to do better and correct for 
toroidal effects as well. Note if the Mirnov coils are mounted just inside a conducting vessel then 
b r = 0, and the second terms are zero. 

Now we must consider a general distribution j(r,0) instead of the line current at r = a. Then using 
complex notation, and assuming a conducting vessel where b r = 0: 



thus relating moments of the external field to moments of the current density. Rogowski coils 
with winding density cos(mG) and sin(mG) would directly measure the LHS of this equation 
(again we knew this from section 9). Now we see a problem: the LHS is an integral over while 
the RHS is an integral over r' and 0. Without some assumptions we cannot work back from the 
externally measured fields to the currents at the surface. This is exactly the same problem that 
we came across in dealing with reconstructing the plasma shape from the multipole moments 
without invoking equilibrium: varying the assumed distribution gives different plasma shapes. 

A typical assumption is 



with jo(r) the equilibrium current, the radius of the resonant surface, 8(r-r mn ) the delta 
function restricting the currents to r = r^. The perturbed field measured at r, produced by the 
current perturbation at r m , is 



Here we have introduced a frequency f = (d/(2n) into the problem: the filaments are rotating. If 
we want to consider b r and be, and work in real space, then 




17.10 




17.11 




17.12 




17.13 




17.14 
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Often we will write an expression for the Fourier coefficients of the field at the resonant surface, 
b mn { r= r mn) = ^ m " r ""y2 ■ Note if we were working on a rectangular vessel we would measure 
for example (see Figure 10.1): 



■j+i 



b„ = b d cos(O)+b rS m(0) = bJ [ -^) cos((m + l)0 + n0-a) mn t) 17.15 



i.e. we see that working with a rectangular contour the phase varies as (m+l)0, instead of m0 
when working with a circular contour. 

Now in toroidal geometry there are complications. We must consider that 
a) jmn produces a stronger field at the inner equator 



b) the perturbation are displaced because of a Shafranov shift 

c) the pitch of the field lines is no longer constant 

We will now derive the toroidal corrections for a case where we are considering surface current 
perturbations for our model equilibrium of section 6. We start with the field line equation for the 
total equilibrium fields, evaluated at a radius r = a p , the plasma minor radius: 

adO Rd(/> 
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B A 



17.16 



With 



R= R g + acos(e) 



1+ — cos(0) 
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Figure 17.2. Field lines around a torus, for the case q = 2. 



We can derive what is sometimes called qmhd, sometimes qi. This is the number of toroidal field 
revolutions a field line must be followed to make one poloidal revolution. Figure 17.2 shows the 
field lines for the case where q = 2. Generally we have 
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Keeping (a/R) 2 terms gives 



2m 2 B. 



00 



r 



MM 



1 + 



( 


Y 


>-( 


a 






u 


J 


V 



2>jl 

2J J 



i + 



( \ 2 f 

a 



1 + 



0.5^ + ^ 



1 



17.22 



The toroidal angle covered when a field line is followed around a poloidal angle of is 
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where we have written q c i rc for the value in a straight cylinder. In this straight cylinder we would 
have the toroidal angle (]) covered in following a field line a given poloidal angle given by 
Equation 17.23 with a/R g = 0, i.e. (]) = q c i rc 0- Now we see that, in transferring to toroidal 
geometry, we must replace the poloidal angle by 0*, where 



123 



Magnetic fields and tokamak plasmas 



Alan Wootton 



= 



IMHD 



0-y(#+^+l)sin(0) 
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i.e. the perturbing fields must have the form (for ~ a) 

(r A m+1 

b ff = b mn [-^j cos{me*+n<t>-co mn t) 



17.25 



Figure 17.3 illustrates the field line trajectories in (<|),0) space. Here we have assumed that 

q MHD (r) = q Q + (q a - qo\~] > w * m 10 the value at r = 0, q a the value at r = a (the plasma edge). 

I 



This allows us to express r = a 1 MHD — — | . In the figure we show examples for a = 0.25 m, 

V q a -q Q 

R g = 1 m, qo = 0.9, q a = 3.2, (3i = 0.5, lj = 0.9. The solid lines are for qMHD = 3.2 (the plasma 
edge), and the broken lines for qMHD = 2. The shear in the q profile is apparent. 




2 Pi 



inside 



outside 

Figure 17.3. The trajectory of field lines in (]), space for q = 2 
(broken lines) and q = 3.2 (solid lines). 



Analysis techniques 



With a set of Mirnov coils spanning a poloidal cross section of a low beta, circular cross section 
tokamak, we can take a Fourier transform in to obtain the amplitude of each component 
cos(m0+n(|)-Q) mn t). If we have a rectangular vessel, then we have shown that the relevant 
expression is cos(m0+0+n(])-G) mn t). This has been done both computationally, and using analog 
multiplexing. We should allow for the toroidal corrections discussed above when performing 
this Fourier analysis; so that is replaced by 0*. Figure 17.4 shows the placement of the coils 
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around a circular vessel. Both poloidal and toroidal arrays are required to determine poloidal (m) 
and toroidal (n) mode numbers. 




Figure 17.4. The placement of Mirnov coils (magnetic pick-up coils) around a vessel 

Figure 17.5 shows an example of data obtained from PBX, a machine with a non circular vessel 
but operated with a circular plasma. Coils 1 through 5 are located on a vertical line at the outer 
equator, with coil 1 the lowest. Coils 6 through 10 are located on a vertical line at the inner 
equator, with coil 6 the highest. Coils 1 1 through 14 are elsewhere; but in particular coil 12 is at 
the outer equator (0 = 0) while coil 13 is at the inner equator (0 = it). Data is shown for a period 
of 10 ms. The scale in volts is shown on the right hand side. A peak first appears on coil 1, and 
then progresses to coil 5. An outside to inside asymmetry (coil 12 to coil 13) is evident. 
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Fz'gwre 17.5. Data from PBX. 



Another technique is to look in the frequency domain. Suppose we have coils placed both 
poloidally and toroidally around the plasma. Then the relationship between the phases of 
different signals identifies m and n without the amplitudes being known. By looking in the 
frequency domain we can reject noise, and other modes at frequencies other than (0^. Taking a 
rectangular vessel as an example (we measure b^), then the fundamental component of the i th coil 
signal will be of the form 



S i = Acos(a) mn t-S i ) 



17.26 



with (flmn identified from power spectra. The relative phase shift 8i will be, for the m,n mode, 

L T r . n _ ,_J 



17.27 



where (Q[,$i) locates the i th coil, 8o+2pki expresses the multi-valued phase property. If we were 
in a circular vessel, then the phase shift would be 
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S: = m 



r 



61 mn 1 



17.28 



By performing a best fit of the measured phase with this expression gives m and n. For example 
plotting Si-n(])i against (to remove the toroidal effects) should give a line of slope m+1. An 

example is shown in Figure 17.6, data from PBX. The best fit to the data was obtained using m 
= 2.n= 1. 




-90 90 180 270 

P0L0IDAL ANGLE 

Figure 17.6. The phase of the experimental data shown in 

Figure 17.5, together with the phase given by equation 17.28. 

with m = 2, n = 1 



Turbulence 



Magnetic coils outside the plasma measure not only the low m,n "Mirnov" oscillations (tearing 
modes), but also higher frequency, higher mode number fluctuations. Many m, n modes are 
possible; Figure 17.7 shows an example where 1 < m < 40, and 1 < n < 12 are considered. 
Each point represents a possible mode combination, the dark area shows those possible when q a 
= 3.2 , qo = 1, and constraints to m are applied (in this case 6 < m < 12). 
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in 




01 2 315 678 9 10 11 12 

Figure 17.7. The (m,n) space, with the limits imposed by qo = 1 
and q a = 3.2. The darker shaded area represents the space 
considered in the analytic model described in the text. 

The measurements of the fluctuating fields are usually restricted to root means square (rms) 
fields, b rrms and be rms , at r > a, the limiter radius. We would like to know how this is related to 

the field components at the resonant surface, so that we can calculate associated island widths. 
Typically we want to determine b rmn from a measured be rms outside r^. There is no unique 
transformation from b0 rms to b rmn , so a model must be invoked. 

We proceed by assuming the fluctuating self generated field be rms measured with magnetic 
probes at r > a is proportional to the required b rmn at ~ a. Evidence for such a model comes 

from correlation measurements between Langmuir probes and magnetic probes; the measured 
berms at r > a is apparently determined by plasma current fluctuations at r ~ a, i.e. at the plasma 
edge. Outside the fluctuating current filaments (r > r^) the magnetic fields can be approximated 
by 



h r =Z( & ™„) 



r 



cos(m^+n^ + ^ mn ) 



17.29 



with § the toroidal angle, the poloidal angle, and a random phase. This ignores toroidal 
effects, which introduce terms proportional to (m+1), and is strictly valid only for small n. If a 
conducting wall is present at r = b then the expression is modified by the factor [l-(r/b) 2m ]. For 
stationary ergodic turbulence the time and spatial averages are the same, so that the root mean 
square is found as a spatial average. Outside the singular surface we measure 
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b rms = b Tjms = b 0iTms = J0.5Xk 



/ \-2(m+l) 

r 



17.30 



To go further and obtain a relationship between be rms at r > a in terms of b rmn at r = r mn we must 
specify 

a given current density or safety factor profile q(r), 
the dependence of b rmn on m and n, 

locate all possible m and n pairs where q = m/n, 
specify upper and lower limits to either m or n. 

In general this must be done computationally. However, some progress can be made analytically. 
We can approximately solve Equation 17.30 for b0 rms at r > a in terms of b rmn at r = r^ if we 
take b rmn independent of m and n over some range to be specified, and a q profile relevant to the 

frV 

plasma edge, for example q = q \ — . Then we can relate the resonance location r mn to a given 

q. 
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With p = r/a defining the measurement position and q = m/n, Equation 17.30 is written as 
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p > 1 is a constant, namely the location of the magnetic pick-up coils used to measure b rms . 

We now replace the summation by integrals. The space we are considering is shown in Figure 
17.4. The first integral, jdn, is from m/q a to m, allowing all n values in between. The second 
integral, jdm, is from mi to m2, with m2 » mi. In Figure 17.4 mi = 6 and m2 = 12. 
Unfortunately our chosen profile allows q < 1 (when r^/a < 1/Vq a ). To restrict n > 1 we should 
choose mi > q a ~ 3. The result of the integration is 
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This expression can be further simplified for p > 1 by noting m2 » 1 (so all terms involving m2 
disappear) and mi > 1 (so the term involving q a ^ ml+1 ^ disappears) to give 



\b rmn \ 2p^ + yq a ln(P) 



17.34 



This is the required result, giving us b rmn at the plasma edge in terms of the measured b rms 
outside the plasma. Because the terms involving m2 have been ignored in deriving Equation 
17.34 it predicts (b rms /lb rmn l)p_ >1 => °°. The more complete expression, remains finite, and 




17.35 



We conclude that the measured rms field outside the plasma (p > 1) is determined only by the 
lowest m number (mi); the higher m numbers fall off too rapidly to contribute. Indeed, if 
Equation 17.34 is taken to represent the results, then a fit of the form b rms oc rY gives the result 

y=m, + l+ — \— 17.36 
21n(p) 

We can also measure the frequency of the fluctuations. We assume that the observed frequency 
is described by an expression of the form 

f ™ ^ A 

17.37 



CO m 



2^ 2m- mn B^ 



\ \e\n dr r ) 



This is consistent with experimental results, with K a factor of order 1. We also know 
experimentally that the term involving E r dominates the term involving K, and that we can 

approximate the expression by 
—cmE 

r - 17.38 

with % a factor « 1.5. E r is approximately constant for p < 1. We now relate r^ to q = m/n by 

assuming a specific q profile. For the simple q profile used in deriving the analytic results the 
transformation from r^ to m and n is given by Equation 17.31. A more realistic q profile gives 
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Substituting either Equation 17.31 (edge model, f = fi) or Equation 17.39 (f = f2) for into 
Equation 17.38 gives an expression for the frequency of a mode of given m and n; 
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These can be normalized to the frequency fm/n=2/l of the m = 2, n = 1 mode, which is often 
measured, to give from Equation 17.40 



J mn 
^ ^ ml fi- 



17.42a 



and from Equation 17.41 
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Figure 17.8 The values of m and n included and excluded by restricting f > 50 kHz, 
when fm/n=2/l = 6.5 kHz. The excluded area is to the left side of each filter line. The 
dashed straight lines represent the limits imposed by 1 < q< 3.2. The dotted curve 
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represents the limits imposed by the frequency filter from Equation 17.42a and 
Equation 17.42b. 

If only modes with f > f^ are considered, we can derive mij m , a function of fmi n /fm/n=2/l an d n, 
which separates modes which are and are not included in any subsequent analysis. The results 
are shown in Figure 17.8 for f^n = 50kHz, fm/ n =2/l = 6.5 kHz. In general, frequency filtering 
should be considered more a restriction on n than on m. 
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18. INTERNAL PLASMA MEASUREMENTS 

So far we have been concerned with measurements of fields taken outside the plasma. In 
comparatively low temperature plasmas (say T e < 50 to 100 ev), we can design pickup coils to 

make internal measurements. There is always the worry that the insertion of such a coil changes 
the very plasma characteristics we would like to determine. This fear is usually allayed by 
monitoring certain characteristic plasma features (sawtooth activity, Mirnov activity, loop 
voltage) to make sure they do not change significantly during probe insertion. Figure 18.1 and 
Figure 18.2 show a possible coil set up which might be used. The coil itself must be protected 
from the plasma, typically by a stainless steel case, possibly surrounded by a carbon shield. The 
geometry of the surrounding materials must be carefully chosen if we are looking at high 
frequencies so as not to cut off the very signals we want to measure. 

Tvpical internal magnetic probe construction 



Single-Coil 
Probe 



Multi- Coil 
Probe 



1cm 



Stainless steel 
stem 



Silica Housing 



Boron Nitride 
Holder 




Sensing Coil cnot to scaiei 

~7U 



Copper Coil 




Boron Nitride 
Former 



Typical probe insertion geometry 
Vacuum seal- 

Flexible 
Bellows 




Figure 18.1. A possible internal coil 
construction. 



Figure 18.2. A typical coil placement, 
showing the flexible belows. 
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Equilibrium 



We first discuss the equilibrium. From the basic field measurements themselves we can, 
assuming circular straight geometry, reconstruct the current from the equations 



1 d 

JU Q r dr 



K) 



jU dr 



18.1 



18.2 



i.e. to obtain j^O") we only need the radial dependence of Be. Unfortunately we have to contend 
with non circularity and toroidicity. One technique which has been applied is illustrated by the 
results shown in Figure 18.3, where small pick-up coils were used to measure the poloidal 
magnetic field at current peak in a small tokamak (TNT inJapan). 



IN- 




OUT 



100G 



The measured poloidal magnetic field at 
the current peak (700j«sec). The dimensions are 
in millimeters. 



Figure 18.3. Equilibrium poloidal fields measured in TNT 

We move to the coordinate system of Figure 1.7. The radial component Br(Ro,z) is measured 
along a vertical line R = Ro, and the vertical component B z (R,0) is measured along a line z = 0. 
The results are fitted to expressions of the form 
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18.3 

n=0 

B z (R,0)=j^b n R" 18.4 



n=0 



The magnetic axis is found where from the zero crossing of the resulting polynomials. The flux 
function is found by integration: 

N 

^(R ,z) = -R Q \Y J ( a nZ")dz + const 18.5 

ra=0 

N 

y/(R ,o)= j^(b n R n+1 )dR + const 18 - 6 

The current density is then obtained as 

MJ = ^-^ 18-7 

" oJ dz dR 

The constants (giving 6Br/5z at z = and 6B Z /9R at R = Rq) must be determined by making 
some assumptions concerning the plasma shape, say that it is mostly elliptic. Some examples of 
the results of this analysis, where N = 5, are shown in Figure 18.4 for the fluxes and Figure 18.5 
for the current density. 
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A rough sketch of magnetic 
ourface at the current maxLnun. 
The equi-psi surfaces are drawn 
for every 7cx 10~ 3 V-sec. Num- 
bers at arrows show the field 
component parallel to the arrows 
in unit of gauss. The shaping 
coil current is presented in 
kAt, and minus sign means that 
the coil current is parallel to 
the plasma current. Coil current 
has up-to-down symmetry. 



Figure 18.4. Flux surfaces reconstructed for TNT. 



BOTTOM 




(cm) 



A typical example of the current density 
profile at the current rising stage (r = 500//sec) and 
at the current peak (f = 700 //sec). Origins of the 
graphs correspond to the center of the vacuum 
vessel. 
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Figure 18.5. Current density profiles for TNT 

Given the current density we can also attempt to deduce the safety factor q. 

Internal magnetic measurements are also used to determine the internal electric field. From 
Faraday's law 




18.8 



18.9 



Figure 18.6. The geometry used to describe the 
measurements of internal electric field 
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Toroidal electric field evolution 



Figure 18.7. The internal electric field prior to a disruption (from Hutchinson) 

Figure 18.7 shows some spatial profiles of E^ from this analysis just before and at a disruption. 
Although the edge electric field goes negative (the negative voltage spike) the internal electric 
field strongly positive. In principle, having measured j and E we could derive the local 
conductivity. 

Mirnov Oscillations 

The same probes used to measure the internal equilibrium properties can be used to look at the 
Mirnov fluctuations (as long as the coils have a sufficiently high frequency response). Data from 
such experiments has isolated the radial dependence of the fluctuating b r , be fields, as shown in 
Figure 18.8. It agrees with our discussion in section 17, namely b oc (rm^r)" 1 " 1 " 1 for r > rmn 
without a vacuum vessel. In the presence of a conducting vessel at r = r w we must account for 
the image currents which flow, so for example we would expect 
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Figure 18.8. The measured radial dependence of the fluctuating poloidal fields 
(Mirnov oscilations) from and m = 2 tearing mode (measurements in TOSCA). 
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19. THE CONDUCTING VACUUM VESSEL 

Skin depths 

There is a vacuum vessel surrounding the plasma. Early experiments used insulating vessels (e.g. 
quartz or ceramic), which were later replaced by metal (stainless steel). In order to introduce the 
toroidal electric field necessary to initiate the discharge, this vessel must be highly resistive. This 
is ensured either by making a thin vessel, using convolutions, or adding at least one insulating 
break. 




. Circuit formation by tfe symmetrical (a) and non-symmetrical 
tb) components of the current in the casing in the presence of transverse 
g»P»- 

Figure 19.1. The currents flowing in a vessel or shell with a poloidal (or 
transverse) gap or insulating break. 




(q) (b) 

Circuit formation by the non-symmetrical component of the 
current in the casing in the presence of longitudinal and transverse gaps: 
fa) longitudinal gaps causing no distortions of the transverse field: 
t» longitudinal gaps distorting the transvene field. 
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Figure 19.2. The currents flowing in a vessel or shell with toroidal 
(i.e. longitudinal) and poloidal (i.e. transverse) gaps. 

There are two types of insulating breaks, as shown in Figure 19.1 and Figure 19.2, (taken from 
Mukhovatov and Shafranov). For transverse gaps (Figure 19.1) any symmetric component of 
current must flow on the vessel surfaces. However, non symmetric components, introduced for 
example by axisymmetric plasma motion, will produce volume currents. At the gap the non 
symmetric currents flow in opposite directions. They will produce a local vertical field, which 
must affect the plasma position, as well as the interpretation of magnetic coil signals. The two 
(surface and bulk) currents have different decay times, as discussed later. 

If both a transverse and a longitudinal gap exists, Figure 19.2 shows that for the non symmetric 
components the placing of the break in poloidal angle is important. A longitudinal gap at the 
outer equator does not distort the field, while a gap at the top or bottom does. 




Figure 19.3. The geometry of a conducting shell or vessel. 

We are concerned with the currents induced in a resistive vacuum vessel lying outside the 
plasma, and their effects on all the measurements we have discussed. As a characteristic 
example, consider a homogeneous field B = Bosin(O)t) parallel to the plane of a conducting plate; 

the inner and outer boundaries of the plate are at a distance b and a from the symmetry plane (see 
Figure 19.3). Maxwell's equations reduce to 



dB 

Vx(VxB) = -// <7— 



19.1 



When the ratio of the plate thickness to the layer thickness is much less than unity, i.e. d/b « 1, 
the penetration of the field is characterized by a time constant 



T = 



jU Q obd 



19.2 



Now we can consider two cases, depending on the ratio of plate to skin depth, d/d^n, where d s kj n 
= (2/((iatt))) 1/2 . If d » d s iQ n , (i.e. at high frequency) we have a thick plate, and we can consider 
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the plate to be perfectly conducting. If d « d s ki n we have a thin plate, and then the characteristic 
penetration time is given by equation 19.2, namely x = (iabd/2. 

The result is the same for the penetration of a longitudinal field into a hollow cylinder, i.e. 
(approximately) a toroidal field into the toroidal vacuum vessel, x is the L/Q time for the decay 
of the poloidal current, with Le = 7t|loa 2 the inductance per unit length, and Q.q = 2na/(od) the 

resistance per unit length. The same result applies for the penetration of a transverse field into a 
cylindrical shell, i.e. a vertical field into a toroidal vacuum vessel. 

When we first initiate some perturbing field, the vessel will first appear as a thick plate, then 
change over to a thin plate after some time x c . This time x a is just the time it takes any dipole 
currents to become homogeneous, and is obtained from x a ~ 7t/(2(fl), where (0 is the value for 
which d s ki n = d/2. We take (d/2) to allowing for penetration from both inside and outside the 
plate (from gaps in the plate), and obtain. 

Ttu^od 1 

t=— 19.3 

16 

Application to a diamagnetic loop 

Here we are considering poloidally flowing currents. Consider an N turn loop wrapped around a 
vessel with minor radius a v , used to measure the plasma diamagnetic current I s = <p±>/B^. After 
a time t > x a we can consider the vessel as thin, and need correct only for the long time inductive 
time constant L/Q.. The voltage around the loop is 

£ dl dl 

— =M ^+L^ 19.4 
N dt dt 

with I s the diamagnetic current, I v the vessel current, M = 7l(loa p 2 and L = 7C(l a v 2 . The vessel 
current is given by 

dl dl 

= n/ v + M^ + L^ 19.5 

dt dt 

where Q. = 27ia v /(od). From 19.4 and 19.5, we have 

r; l i 



/=— \edt + — e 19.6 

nm[{ a 

Therefore we can correct the measured voltage £ to obtain the required diamagnetic current I s 
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Application to position measurement 

Here we are considering toroidally flowing currents. If we are using the moment coil method to 
measure the plasma position, then the coils measure the current center of any currents within the 
contour on which the modified Rogowski and saddle coils are wound. Therefore, if these coils 
are outside the vacuum vessel, they will be sensitive to any currents induced in the vacuum vessel 
itself. If we consider times t > x s , then we need only worry about the homogeneous vessel 
currents, driven by any non symmetric part of the flux function, *P a . A vessel current with 
equivalent surface current density 

ad dm 

i=oEA=^ L ^ -i 19.7 

v v * v 27lR v dt 

will appear, where subscript v refers to the vessel. This asymmetric part of the flux function can 
be taken, for example, from the analytic expressions given in section 6. 

Another method of calculating the induced vessel (toroidal) currents is to represent the vessel as a 
number of filaments, as discussed in section 1. Induced currents are then calculated through a 
mutual inductance matrix. The plasma might be specified analytically or as a single, or number 
of, filaments. Each vessel filament must be given a minor radius, which might be taken as d v /2, 

half the vessel thickness, so that the self inductance is finite. Even the non homogeneous current 
components can be calculated in this way, if enough filaments suitably spaced and connected are 
included. 
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20 . THE IRON CORE 

Many tokamaks have an iron core to ensure good coupling between the primary winding and the 
plasma. This iron core has the additional advantage of keeping 'stray' fields away from 
diagnostic apparatus. However, it complicates the study of the plasma equilibrium, because the 
free space expressions for the fields produced by circular conductors are no longer applicable. 
Here all I want to do is to point out that iron is a pain to deal with. 




Figure 20.1. A straight iron cylinder Figure 20.2. An iron core with an air 
surrounding a plasma gap, with a linked plasma (current I). 



The iron core introduces a new boundary condition. It is often stated incorrectly that the lines of 
force enter a medium of infinite permeability (the iron) perpendicularly, but this is not always 
true. A simple example of a straight iron cylinder surrounding a straight wire with current shows 
this not to always be the case: in Figure 20.1 the field lines are tangential to the iron boundary. 
A correct boundary condition, for the normal component of the induction, is 

B n (air) = ju Q H n (air) = juH n (iron) = B n (iron) 20. 1 

For permeability u very large we can take H n approaching just inside the iron, so that 

H n (iron) = 20.2 
(from the continuity of B n ) and add the boundary condition 

H T is continuous 20.3 
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In reality (I is not infinite, but the limit works well. 

First consider an air gap inductor, as shown in Figure 20.2. We assume that inside the iron there 
are laminations which ensure no current, so that 

V x H = 20.4 

Then H is derived from a single valued potential given by Laplace's equation 

VV = 20.5 

and the boundary condition on B n (Equation 20.2) is equivalent to 
aw 

= in the iron. 20.6 

an 

The only possible solution of equation 20.5 and equation 20.6 is \|/ = constant . Therefore H (but 
not B) must be zero inside the iron. Continuity of the tangential component of H then shows that 
the lines of force in the air must be perpendicular to the iron. Therefore in air the magnetic fields 
are given by the usual equations, together with the boundary condition 

e„xH=0 on the iron surface. 20.7 
In the iron H = 0, but B is finite. Since VxB = in the iron, B(iron) = V$, and $ satisfies 



V 2 4> = 

-r-=B n (air) 
an 



20.8 
20.9 



which is known. 



contour 1 





iron 




y 


> i « 
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Figure 20.3. The geometry of an iron core and 
a linked plasma of current I. The contour 1 
considered in the text is shown. 

Now consider the case shown in Figure 20.3, an iron core with a central limb as on TEXT. If the 
total number of turns linking the center core is 0, and there is no current in the iron (Equation 
20.4) then jH-dl = for any contour 1 inside the iron. Then H derives still from a single 
potential, and all our previous steps are valid. However, if a finite number of amp turns I links 
the iron, then jH-dl = I. This means that we can still write the field in the iron as V4>, but now 4> 
is multi valued, increasing by +1 once around the contour 1. That is, in this case, Equations 20.8 
and 20.9 apply in the iron, but with the additional constraint 

£v4>«dl = / 20.10 

In this case, had we assumed perpendicularity, we would have obtained jH-dl = 0, while in reality 
jH-dl = I. Thus the field lines only enter the iron core perpendicularly if no current flows in the 
iron, and the net ampere turns is zero. 

The only place that the iron really affects magnetic diagnostics is in the equilibrium 
reconstruction. No longer can we use the free space expressions for a circular current filament, 
but they must be modified to satisfy the boundary condition. In toroidal plasma devices we 
almost always satisfy the conditions necessary for the boundary condition Equation 20.7. One 
way to model the effect of iron is by placing additional circular filaments inside the iron itself, 
with currents chosen to satisfy Equation 20.7 at a given number of locations. This is illustrated 
in Figure 20.4: the 'image" filaments must have a current flowing in the same direction as the 
filament in air, so that there is an attraction between the filament in air and the iron itself. As the 
filament in air gets closer to the iron, the "image" current must increase, and so the attraction 
must increase. This means that an iron core can lead to an axisymmetric, or n = 0, instability. 
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z 




www \Wi iron boundary 

where exH n = 



filament in air 



+ R 



Figure 20.4. An iron core and plasma current (filament), together 
with the image currents necessary to ensure the boundary condition 
are satisfied. 



In principle the method of images can model a true three dimensional iron core, although it 
would not help much because most of our equilibrium work involves two dimensions only. 
Therefore we often assume a two dimensional iron core which is produced by the toroidal 
revolution of the actual core. In this case there are analytic expressions for the additional field 
components produced by the iron. 

Because the iron core is really three dimensional, it introduces a perturbation with a toroidal 
mode number n equal to the number of outside return limbs, Ni. As such any fields, such as the 
vertical maintaining field, will have an n = Nj component. Calculations show that this 
perturbation can be as much as 10%, which might be thought to introduce a toroidally dependent 
plasma position. It does, but only by a very small amount. Suppose the toroidally symmetric 
vertical field B z has an additional perturbing component b z cos(Ni(|)). From the field line equation 

this variable field component will displace lines of force according to 



dz Rd(f) 



20.11 



b z cos(N^) B t 



or 




20.12 
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Therefore the maximum displacement is 

A^max = -— 2 - 20.13 



which is typically a few mm. 
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21. TOKAMAK POSITION CONTROL 

The axisymmetric instability 

We first calculate the growth rate of an axisymmetric instability in a tokamak. We only consider 
vertical motion, because it is easier than the calculation for horizontal motion. The driving force 
is written in terms of the decay index n = -(R/B z )(dB z /dR). We consider a tokamak surrounded 

by a conducting vacuum vessel. There is assumed to be a transverse (poloidal) insulating gap in 
this vessel, so that only dipole currents can flow. The equation determining the vacuum vessel 
current I s is 

— (LI )+— (M I n 21.1 

dt y 5 s dt y sp p} s 5=0 

Here I p is the plasma current, L s is the vessel inductance, Q, s the vessel resistance, M sp the 
mutual inductance between plasma and vessel. We introduce the vessel time constant x s = L s /Q s . 
We can approximate the vessel as a circular shell, so that 

L 21.2 



4 



where R s , r s , 8 S are the vessel major radius, minor radius, thickness, and a is the conductivity. 
Approximating the plasma as a filament initially centered within the vessel (R = Ro = R s , z = z s = 
0), the mutual inductance and its spatial derivative are given by 

M 21 .4 

P 2r 



dz 2r 



21.5 



From equations 22.1, 22.2 and 2.5 we can derive the relationship between the dipole current in 
the vessel and the plasma displacement z: 

^ + i = - 2 i* 21.6 
dt -r 7Vr„ dt 
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The equation determining the plasma displacement z is 



m- 



df 



~0 = -2nR p I p [B R +B s + AB] 



21.7 



with Br the equilibrium major radial field, B s the major radial field from the vessel currents, and 
AB a perturbation to Br. We can take the mass m of the plasma to be zero in the presence of the 
vessel; we will find the vessel slows the motion sufficiently for the m(d 2 z/dt 2 ) term to be 
vanishingly small. Using 



R R 



B = 



4xR 



21.8 



21.9 



T = ln 



v r p J 



2 Hl 2 



5 I 4r 



o 



21.10 



21.11 



where r p is the plasma minor radius, and we assume R s = R p = R. Taking AB = 0, equation 22.7 
can be written 



= - 



AnR 1 



+ • 



4r 



21.12 



Equations 22.12 and 22.6 define the problem; they have solutions 



z=z e 



I =1 e* 



21.13 
21.14 



The growth rate is 
1 n 



7 = 



n. = 



T s n + n s 
2R 2 



k r 



21.15 



21.16 
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From equation 22.15 we see the plasma is unstable if n < 0. 

Consideration of in-out motion is more complicated, because we must conserve poloidal 
magnetic flux. This is ensured by using the equation 



-[L p I p + *R 2 p B z+ M ps I s ]=0 



21.17 



where the plasma inductance is 



In 



v 



V r P J 



■2+^ 
2 



21.18 



The unstable motion we nave derived for a plasma with decay index n < must be controlled by 
feedback. Figure 21.1 shows such a system. Sensing coils might be the multipole moment coils 
discussed in an earlier chapter, or single coils. Sensors can be placed inside or outside the vessel. 

^ _ _ -(^) feedback winding 




feedback winding 
Figure 21.1. Geometry of a feedback system 
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Analysis of sensor coils allowing for vessel currents 

If single coils were used, the analysis used to determine plasma displacement z is to take the 
difference signal, and divide by the sum signal. For a positive plasma current (into the plane) 
and a positive vertical plasma displacement z, the upper coil signal is Bgi = (j,oIp/(27l(r c -z)) and 
the lower signal is B02 = (JoIp/(27t(r c +z)), so that for z « r c : 



B m +B d2 



1 ^ 



r c - z r c + z 

1 1 
+ 



\r-z r+zj 



— 21.19 



We now consider what happens when we allow for the presence of vacuum vessel currents I s and 
feedback currents If, both of which produce fields seen by the sensing coils. 

The dipole model 

The feedback windings are represented by a shell of radius rf and thickness 8f with a dipole 
current distribution, with current (area) density j^f = j ( | ) fosin(0). Then 

h = \ i« r f s f d0= lr f S fko 21 - 20 



If is positive when the upper winding current is into the plane. This feedback circuit produces a 
major radial field for vertical position control, and (r c < rf): 

B t = K r I f , Kf = — 21.21 

J J J J A r 

<+r f 

That is, a positive If produces a signal Bqi f = Bf (which will be negative for positive If) in the 
upper sensing coil and a signal Be2,f = -Bf (which will be positive for positive If) in the lower 
sensing coil. 

The vacuum vessel is also assumed to be a shell, of radius r s and thickness 8 S , carrying a dipole 
current j^ )s =j(j )s osin(0). Then 

I s =\j, s r s S s de = 2r s Sj, s() 21.22 

o 

This produces a major radial field for sensors inside the vessel (r c < r s ) 
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B SJ = KJ S , Ks^-f 1 21.23 

Therefore the signals seen by the sensing coil are Bq\ s j = B s j in the upper coil and Be2, s ,i = -B 
in the lower coil. If the sensing coils are outside the vacuum vessel (r c > r s ), then the maj 
radial field is 

B s , e = K s , e I s , K s,e=7T 21.24 

Therefore the upper coil sees a signal Bq\ s e = B s j, and the lower sensing coil a signal Be2, s , e = 
B s e . These results are summarized in table 1. 



Table 1. Fields seen by the Be sensing coils, 
dipole current model. 



plasma 
feedback 

vessel 

inside 

plasma 

feedback 

vessel 





B 82 


(r c >r s ) 






Moh 


2n(r c - z) 


ln(r c + z) 


Mo 1 / 


Mo 1 / 




Ar 


Mo 1 /* 


Mo 1 /, 






(r>r c ) 




Moh 


Moh 


2n{r c - z) 


ln(r c + z) 


Mo[ f 




Ar 
<*r f 


Ar f 


Moh 


Moh 







The apparent displacement z app is derived using the equation 

5^ = |X-|X = z_ +a L + a k 2L25 

r e E*« + E**2 r c s I p f I p 
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where the values of a s and (Xf, derived from the expressions in table 1, are given in table 2. 
Without the correction factors a s and 0Cf, z app = z. Moving the sensors from inside to outside 
changes the sign of cc s . When the plasma (I p > 0, into plane) moves upwards (z > 0) it will 
induce a negative I s . For sensors outside the vessel, the vessel currents then produce an apparent 
displacement smaller than the real displacement (i.e the vacuum vessel shields the sensing coils). 
For sensors inside the vessel, the apparent displacement is larger than the true displacement. The 
feedback windings complicate this simple process. If we wanted to, the effects of the feedback 
coupling to the sensors (i.e. finite 0Cf) can be removed by adding to each signal a term 
proportional to If itself. This cannot be done for the coupling from the vessel to the sensors, 
because we do not measure the vessel dipole currents. 
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Table 2. Correction coefficients cc s and (Xf 
for Be sensors. 



correction 
factor 



outside (r > r ) 



7Tr 
2r 



a, -^<0 



2r f 



inside (r s > r c ) 



For comparison, table 2 shows the coefficients oc s and 0Cf for the multipole moment coil system 

[1], with coils either inside or outside the vessel. Note that the multipole moment sensing coils 
are only sensitive to currents inside the contour on which they are wound, and thus (Xf = 0. 

Table 3. Correction coefficients oc s and (Xf 
for multipole moment sensors. 

correction 
factor 

outside (r c > r s ) 

7lr 

a s -->0 

a f 
inside (r s > r c ) 
a s 
a f 

Numerical values applicable to TEXT are given in table 4, for r c = 30 cm (inside) or 34 cm 
(outside), r s = 32 cm, rf = 50 cm. We see that, as far as pickup from the shell currents (a s ) is 
concerned, the Be coils outside the vessel are no worse than the moment coils. 
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Table 4. Numerical values of the correction factors Oj and (Xf. 

R = 100 cm, r c = 30 (inside) or 34 (outside) cm, r s = 32 cm, rf = 50 cm 



dipole model, multipole coils outside 
dipole model, multipole coils inside 
dipole model, Be coils outside 
dipole model, Be coils inside 



-1.67 







1.48 



1.48 



-0.94 







-1.07 







The feedback model 

Here I want to illustrate how a feedback system is analyzed. We take the mass-less plasma 
model, in which the plasma is always in equilibrium due to currents in the vacuum vessel and the 
feedback windings, and the (unstable) driving term Br. The vacuum vessel (i.e. conducting 

shell), feedback windings and plasma all having the same major radius R. The equations 
describing the plasma motion z, shell current I s , feedback current If and feedback voltage Vf are : 



That is, in open loop the plasma displacement induces a dipole shell current as well as a dipole 
feedback current. The dipole currents in the shell and feedback circuits are themselves coupled. 
In closed loop a proportional (g) and derivative (gtd) gain term are considered The mutual 
inductances My and self inductances Lj between the various circuits (i,j = p for plasma, s for 
shell, i.e. vacuum vessel, and f for feedback) are: 






21.26 



Z=z Q -z i 



■app 
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M, 



sf 



4r f 



* 2r f 
M 

4 



21.27 



4 



In — +--2 



a p 2 



g is the linear gain in the feedback circuit, t a is the time constant of the servo system, td is the 
time constant of the derivative gain, t2j is the resistance of the circuit i, zq is the required 
position. The equilibrium major radial field Br is given in terms of the equilibrium vertical field 
(-B z ) and the decay index n =-(R/B z )(9B z /6R) as 

B =^ 



R 

z 4xR 



21.28 



T = ln 



f \ 
SR 

\ a nJ 



+ -+ Pt-- 
2 Hl 2 



For vertical (up/down) instability without feedback control n < 0; following a vertical 
displacement the plasma then experiences a major radial field Br in a direction to cause further 

plasma displacement. We can model any unstable motion through the choice of n: without 
feedback but with vacuum vessel n should be chosen such that the measured growth rate is: 

1 



ft, = 2- 



n 

t s n + n s 
R 2 



21.29 



r 2 r 



Note that the vertical field required for major radial equilibrium is -B z as defined above. We 
normalize the equations, letting 
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21.30 



Values of x s and Xf are derived from the values of oc s and ccf in table 4. Representative values for 
TEXT are given in table 5. 

TABLE 5. values of x s and xf used in the modeling 

x s x f 
multipole outside 0.5 
multipole inside 
Be outside 0.5 -0.3 

B e inside -0.5 -0.3 



The normalized equations are: 

-n£+ k s i s + k f i f = 
d£ di . di f 



dC 



di, di f i 



jU f — + m fs —^+—^+- L = ^- 
dt ' dt dt T f T f 



dv f g 

f a dt k, 



21.31 



(fo~C- xj s - xj f )- T c-^(Co -C- x s i s - xfy ) 



These can be solved either numerically or by Laplace transforming. These equations are then 
Laplace transformed (a superscript "~" denotes a Laplace transform, s is the Laplace variable) to 
give 
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-nC(s) + kj s (s) + k f i f (s) = 

jU s s£(s)+(s + 1X0) + m sf si f (s) = 

~ ~ I 1 ] ~ v f (s) 21 32 
jU f s^(s) + sm f JXs) + \ s + — / f 0)-- i — = " 

v f (s) = j- | (j) -(1- T z s)^(s)+ x,i,(s)+ */}(*))] 

where X£ = x a + Xd- Solutions for the displacement ^, vessel current i s , feedback current if and 
feedback voltage Vf can be obtained; for the displacement C, is given by: 

C(s) A+Bs 



Q(s) D+Es+Fs 2 



21.33 



where 

A = gk f 

B=-g(k s m sf -k f )=0 
D = k f (n+ g)+ ngx f 

E= -gx s m sf {n + k s ji s ) + k f ^i + k s /x s + (n + A: / // / )r / ]- gk f t L + gx f (n + k s m s - nt s ) 
F= m^x,T z + T f k s (l- m sf m fs + kju s ]- t s x f g(n + k s m s ) 

21.34 

The term appearing due to the coupling between sensors and feedback windings are in bold. The 
remaining terms are identical to those previously derived, except that the values of x s to be 

considered are different (negative values are now permitted). 

For stability, the Routh criteria can be applied, which requires 1) D > 0, 2) E > 0, and 3) F > 0. 
These three criteria define the available operating space in (g,X£) space for a given instability (i.e. 
for a given decay index n). 

Application of stability criteria. 

Values relevant to TEXT Upgrade are t s = 7ms, tf = 137ms, Xf = 20, m s f = mf s = 0.67, k s = 2.26, 
kf = 1.5, (l s = 1.59, (if = 1.06, and the pairs of values for x s and Xf given in table 5. Experimental 
observations of open loop vertical motion show Xi ns tability = 1/Y ~ 100 ms, so that a suitable 
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value for the decay index in the modeling is n = -0.2. We now compute the stable regions in 
(g,T£ space), the results of which are shown in Figures 21.2 through 21.5. For the values 

considered the best system is that of multipole sensing coils inside the vessel (Figure 21.2), 
closely followed by the single Bq coils inside the vessel (Figure 21.4). For sensing coils outside 
the vessel the operational space is restricted, more so in the case of single Be sensors than for 
multipole sensors. This is seen by comparing Figures 21.3 and 21.5. 

50 
40 
30 
20 
10 




-3-2-1 12 3 -3-2-1 12 3 

Figure 22.2. Moment coils, inside vessel. Figure 22.3. Moment coils outside vessel. 
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Figure 22.4. Single Be coils inside vessel. Figure 22.5. Single Bq coils outside vessel. 
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22. MAGNETIC ISLANDS 

We consider the effect of a helical perturbation on a tokamak equilibrium. We will not consider 
anything other than vacuum fields, i.e. we will not allow the plasma to respond to the applied 
fields in any way. I have adapted the work of S. Matsuda and M. Yoshikawa, in Japanese 
Journal of Applied Phys. 14 (1975) 87. 

The field line equation is 
dl 

— = const. 22.1 
B 

i.e. working in a quasi cylindrical system based on the plasma center 
dd RB e 1 

— = - = =i(r) 22.2 

d<t> r q(r) 

q is the safety factor, i the rotational transform. Expanding i(r) near r = rg, the radius of the 
resonant surface, we obtain 

dd / dl{r) s \ < ^ _ 

— ~l(r Q ) + -^(r-r Q ) = l Q + lx 22.3 
d<p dr 



where x = r - ro, io = i (ro), and i = 

dr 

have 



. To lowest order in the perturbing field we also 

r=r 



rd6 rdO B e (r) 

where b r (r,0,4>) is the radial component of the error field. We ignore the component because it 
is small compared to the equilibrium field Be. Assuming small islands, we also ignore any radial 
dependencies of b r and Be. b r is expanded as a Fourier series 



a mn s>in(m0 - n0) + a m , n s'm(m0+ n<f>) 
+c mn cos(mO- n0) + c' m ,„ cos(m0 + n0) 



22.5 



Because we shall see that only resonant components matter near r = vq, we keep only the first and 
third terms which satisfy 
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m 1 
n l n 



Then 



b r (r ,e,(f) = b mn sm(m0 -n</> + 0) 



Next we define a new variable 
rj= mO-n(j) + (3 

so that equation 22.4 becomes 

dx _ b mn sm{ri) 
r Q d0 B d (r Q ) 

We also obtain from equations 22.3 and 22.8 
drj n 



dO 



+ m 



l + ix 



Substituting equation 22.9 into equation 2.10 gives us a non linear equation 

r l b sin(7) 



d T] _ d_ 

dO 2 ~~ n do 



1 


r 






l- 








— X 




Jo 


V 


h J 





= n- 



B d (r Q ) 



= -Asin(/7) 



where 



A = —n 



l l B d( r o) 



22.6 



22.7 



22.8 



22.9 



22.10 



22.11 



22.12 



We can take A > for convenience. We see the behavior of the field lines is analogous to the 
motion of a particle in a periodic potential U = A (l-cos(T|)), when we regard T] and as position 
and time. Multiplying equation 22. 1 1 by dr]/d0 and integrating over gives 



22.13 
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with E- 



ifdnV 2 E 

—I — a constant of integration analogous to a kinetic energy. Let k = — , so 
2\dQJ 2A 



that 



-de) 



k 2 - sin 2 ! — 



22.14 



Two types of solution are possible, as shown in Figure 22.1. If E < 2A the field lines move 
periodically within a limit of r|. Field lines with E > 2A are "passing". 




Figure 22.1. Behavior of field lines in phase space. 
The maximum excursion is found when E = 2A, or k = 1. Then from equation 22.9 we have 



w 



= \dx = \ 



B(r \dH 

D ayi 



22.15 



dO 



Using equation 22.14 the integral in equation 22.15 is found to be -j= , so that we finally obtain 



w =4r n 



I 



1 K,„ 


I 


mB (r o ) 


-r Q l 



= 4 



22.16 



dr 
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23. SOME EXPERIMENTAL TECHNIQUES 

Coils winding 

The rogowski coil is simply made by obtaining a delay cable, and returning the wire down the 
center of the delay line (to ensure no net single turn is left). More complicated coils must be 
made by using variable winding densities (i.e. changing the pitch) or varying the cross sectional 
area of the former on which the coil is wound. 

Interference suppression 

Electrical equipment designed to produced RF energy such as generators, and switching 
phenomena in electrical circuits, create RF spectra which must be contended with. 
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Figure 23.1. An illustration of interference paths and suppression techniques. 

The sources of interference are illustrated in Figure 23.1. The interference propagates either 
down lines (cables) or by direct radiation. If the wavelength is large compared to the dimensions 
of the interference source only minor radiation will result, which is mostly found along the lines. 
This is the case for frequencies up to 30 MHz. When the dimension of the interference source is 
about that of the wavelength the interference energy will travel by radiation. The dominant 
frequencies are those where the interference source are 1/4 or multiples of it. Favorable radiation 
conditions imply reduced line propagation (because of increased line attenuation) . Therefore the 
two propagation paths, comprising direct and capacitive or inductive coupling, suggest two 
means of suppression, either line attenuation or de-coupling attenuation. Line attenuation is 
effected by filters. De-coupling attenuation is effected by the construction of the sensor coil and 
the associated connecting lines. 

A common problem with probes is capacitive pick-up. To test for this pick up on simple sensor 
coils, two identical and adjacent coils can be connected in series. Depending on the orientation, 
the signal obtained should be twice that measured with a single coil, or zero. If the coils 
connected in opposition do not give a zero signal, then capacitive coupling effects should be 
considered as a possible source of error. Capacitive coupling can be over using a grounded 
screen or can around the sensor. 

Screened rooms 

The requirement is to screen a room in which a sensitive measurement is being performed from 
external interference, or to accommodate apparatus which radiate interference in a screened room 
to keep the surroundings free from interference. The basic method is to use cages of wire mesh, 
or metal sheet. Both electric and magnetic field components must be considered. Units used for 
effectiveness are the decibel : 
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The wire mesh works to screen electric fields because the external flux lines mainly end on the 
mesh. The effectiveness depends primarily on the size and type of the mesh. Magnetic screening 
is effected by induced currents; DC magnetic fields are not screened, and low frequency AC 
magnetic fields are only poorly screened by non magnetic materials. With increasing frequency 
the magnetic shielding improves and approaches a finite value. Double screens, insulated from 
each other except at one point, improve the screening. These rooms work well to 20 MHz. 
Above this the screen room size can equal the cage dimension, causing resonances. 

Sheet metal rooms have better screening properties than double walled wire mesh, but breathing 
is a problem. The screening against electric fields is ideal since no flux can penetrate. The 
screening of the magnetic component improves with increasing frequency due to the skin effect. 

Honeycomb inserts are also used. The grids are wave guides (with the frequencies considered 
below cut-off), the screening effectiveness of which depends on the ratio of depth to width of the 
honeycomb up to cm wavelengths. They are used for 100 kHz < f < 1000 MHz. 

Misaligned sensor coils 

Typical tokamak requirements include the measurement of poloidal fields in the presence of a 
much larger toroidal field. A small misalignment of the coil will then introduce unwanted field 
components. There are a number of solution to overcoming this problem 

a) subtract data obtained with only the (unwanted) field component by energizing only the 
offending windings 

b) make use of the differential nature of a pick-up coil signal. For example, consider the 

toroidal field to be the offending field, so that the pick-up coil measures 
dB = dB poloidal + dB toroidal ff ^ toroidal field is almQst steady state (d/dt ^ Q) during the tim£s Qf 
dt dt dt 

interest, then the differential signal during this time is approximately that required (i.e. from the 
poloidal field component only). Therefore the temporal integration should be started as late as 
possible. 
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